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ABSTRACT 



A single period inventory aodel is considered for 
supporting an overhaul depot. Demand for any stocked item 
is assumed to come from several repair lines, each of which 
has a different demand distribution and unit shortage cost. 
An approximation for the expected shortage costs was 
proposed by McHasters. This thesis evaluates the accuracy 
of that approximation as well as two additional ones 
proposed by the author. For the case of independent Poisson 
distributions two of the approximations are shown to give 
the exact expected shortage costs. Simulations were needed 
to evaluate the accuracy of the approximation formulas for 
binomially and normally distributed customer demands. By 
imposing restrictions on the range of certain parameters, 
acceptable levels of accuracy can be achieved for the two 
approximations. The imposed restrictions leave sufficient 
freedom to then determine the optimal quantity to stock for 
a given item. 
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I. INTRODUCTION 



A. GENERAL CHAR ACTERISTICS OF INVENTOR! MODELS 

Inventory models have been developed and used to deter- 
mine the optimal policy for running inventory systems since 
1915. The objective is usually to minimize the annual 
expected variable costs associated with running the inven- 
tory system. Usually the costs are aggregated under three 
cost categories: 

1. The order costs depend on the frequency of orders for 
resupply. 

2. The holding, cos ts are those costs associated with 
carrying stock in the inventory system. 

3. The shortage cos ts are those costs which are incurred 
whenever the system runs out of stock. 

When a specific model is constructed and applied to a 
real world situation two additional goals tend to be sought. 
The first objective is to reflect within the model the rele- 
vant features of the specific type of real world situation 
to a satisfying degree. The second objective is to limit 
the costs associated with applying the inventory model. 
These costs consist mainly of costs of data collection, data 
transfer, and computation. In general, a better result in 
attaining one goal implies a reduction in fulfilment of the 
other goal. 

This thesis will consider that special class of inven- 
tory models known as the single period model. Such a model 
attempts to determine an optimal quantity of an item to buy 
and stock which will minimize the holding and shortage costs 
over the period’s duration. Order costs are usually ignored 
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iecause an order is assumed to be placed at the start of 
every given period. 

B. DETERMINING THE SHORTAGE COST 

Of the three cost categories the shortage costs are 
usually the most difficult to express in dollar values 
because the data needed to evaluate such costs is seldom 
available in a readily useable form. It typically is 
included as part of overhead costs or indirect costs. 

Assuming that shortage costs can be quantified, the 
usual approach is to associate the event of ’no delivery' 
with a fixed number of dollars for each item not delivered 
immediately upon request. This can be refined by adding a 
time-dependent dollar cost which is a linear function of the 
length of the time period between the occurrence of the 
shortage and the end of the period. An implicit assumption 
is that the customer population is homogeneous. 

This thesis, considers the case where the customer popu- 
lation is not homogeneous. In this case several different 
major customers are assumed and there is a different demand 
distribution and fixed unit shortage cost for each major 
customer. The unit shortage cost is assumed to be constant 
for each customer. 

C. THE NARF MODEL 

The specific model considered in this thesis is the 
single period model presented by McMasters [Ref. 1]. That 
model is intended to provide supply support for overhauls of 
repairable aircraft systems done by Haval Air Rework 
Facilities (NARF’s). The following assumptions apply: 

1. Within the modelled single period, repair parts are 
needed for the overhaul of different end items. Each 
end item can be viewed as a different major customer. 
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a demand for a 



During the repair of an end item, 
given part occurs with probability p, where p is 
greater than zero , ( and no demand will occur with 
probability (1-p) ). Demands are satisfied on a 

"first-come, first-served" basis. 

2. The repair cycles for each type of end item are 

nonoverlapping and one cycle immediately follows its 
predecessor. There may be a different number of 
repair cycles per period for different end items. 

3. A fixed unit shortage cost is associated with a 

specific end item. This unit shortage cost may be 
different for different end items. 

4. The quantity to buy for inventory is computed and 

placed in stock at the beginning of each guarter 

(period) and no replenishment during tne guarter is 
allowed. It is assumed that the supporting Naval 

Supply Center (NSC) maintains the stock or will 

obtain it by the beginning of the guarter when it is 
needed. 

Within this setting the expected shortage costs for the 
period are a function only of: 

1. The customer specific fixed unit shortage cost; 

2. The order in which demands occur during the period; 

3. The initial stock available to satisfy demand. 

When a demand occurs, the amount demanded is one unit 
and the order in which each demand occurs depends on tne 
probability distribution of the inter-demand times for each 
customer, where the inter-demand time is the time between 
successive occurrences of demand by that customer. The 
total demand by each customer during the period is binomi- 
ally distributed and is assumed to be stochastically inde- 
pendent of the demand of any other customer. 
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McMasters suggested, an approximation formula for the 
expected shortage cost per period (two customers) [Ref. 1], 



expected shortage 
cost per period 



where 
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The first term attempts to approximate the average unit 
shortage cost. The second term evaluates the expected 
number of shortages per period. The first term represents a 
convex combination of and 77^ weighted by ^ and //> . Its 
range of values is bounded by /n, and 7^. A large value of 
drives the value of the first term towards whereas a 

large value of drives it towards 77 ,. 

Several theses have expanded the work of McMasters. 
Slaybaugh [Ref. 2] analyzed the demand data and the rework 
schedule during 1980 for the TF-34 aircraft engine being 
reworked by RARF at the Naval Air Station Alameda, 
California. Combining those data with the engine’s bill of 
materials, he was able to estimate the probability of 
replacement per repair action and conducted a parametric 
analyses to study the relationship between the shortage and 
surplus costs. 
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Asselin [Ref. 3] applied dynamic programming in order to 
expand the single period model to a maiti-period model. The 
development of a "good" formula for the expected shortage 
costs per period in the single period model is a basic 
building block for an efficient multi-period model. 

D. PURPOSE 

The purpose of this thesis is to evaluate the approxima- 
tion for the shortage cost term suggested by McMasters and 
to investigate alternatives to that approximation. 

E. PREVIEW OF ANAL! SIS 

Chapter II shows that Monasters' approximation formula 
provides the exact expected shortage costs per period in 
case of an arbitrary number of mutually independent Poisson 
distributed customer demands. Because Monasters' approxima- 
tion formula does not provide the exact expected shortage 
costs per period in case of binomial distributed customer 
demands, other approximation formulas are presented in 
Chapter III. 

Simulation is then used in Chapter IV to VI to evaluate 
the performance of the different approximation formulas. 
Chapter IV describes the simulation methodology. Chapter V 
contains simulation results and discussion for the case of 
two customers each with binomially distributed demands. 
Chapter VI expands the analysis to more than two customers. 
The normal distribution is used to approximate tne sum of 
the binomial demand distributions and McMasters’ approxima- 
tion formula is applied. This procedure is then analyzed by 
running simulations. 

A final summary with conclusions and recommendations is 
presented in Chapter VII. 
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II. POISSON DISTRIBUTED SOSTOMER DEMANDS 

A. PRELIMINARY REMARKS 

Suppose that demands for each customer can be described 
by independent Poisson processes. In general, the Poisson 
distribution is often used to model infrequent random 
events. For instance, the UICP inventory model of the 
US Navy for consumable items assumes the very low demand 
items to be Poisson distributed ( an average demand of 1 
unit or less per year ) . 

If the Poisson assumption is valid then the expected 
shortage costs per period are exactly equal to HcMasters' 
approximation and can be easily derived using conditional 
probability arguments. Because McMasters' approximation 
formula basically multiplies an average of the unit shortage 
costs by the expected number of shortages per period it will 
be named ' monovariable approximation formula’ (i.e. only the 
sum varies) . 

B. COMPUTATION OF CONDITIONAL EXPECTED VALUE 

Ross [Ref. 4 pp. 75-76 ] determines the conditional 
expected value of given X t * ... + = x, wnere 

Xj _, ... fX^ are independent Poisson random variables with 

means , t/ U n • 1 

First the conditional probability that X 1 = k given 
Xj_ + ... + X n = x is calculated: 



iRoss computes the conditional expected value for the 
case of two Poisson distributions. His path of reasoning is 
applied here to compute the conditional expected value in 
case of n Poisson distributions for any n > 2. 
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Therefore the conditional distribution of X L given 
... + X n = x is the binomial distribution with pararae- 

. By symmetry it follows that rhe 
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conditional distribution of X i given X x + ... ♦ X^=x is the 

binomial distribution with parameters x and — \ ^ 

for i = 1, ... , n. The associated conditional expected 



values are : 



E[ X' | X x + . . .+X 0 =x ] = x 
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C. DERIVATION OF THE FORMULA FOR EXPECTED SHORTAGE COSTS 
PER PERIOD 

Given a shortage occurs, the demand generating the 

//■/ 

stockout is due to customer i with probability — — - and the 
penalty is 77/ dollars. The expected penalty cost for the 
stockout is therefore 

77 /—. dLL 

1=1 /' x + **• 

Now the number of stockouts given x demands in the period is 

'o if x<Y 

x-Y if x> Y . 

The expected number of stockouts is therefore, 
on 

^ (x-Y) P (x) , 

x= Y+ 1 

where P (X) denotes the probability mass function of X. 

Each of these stockouts results in an expected cost of 



_ _ /<V_ 

jMi * 

Therefore tne total expected shortage costs per period can 
be written as: 



03 

v~ 



/<• 



(x-Y) P(x) 



x=Y + 1 

The last expression is exactly the discrete analog to 
the approximation formula derived by McMasters [Ref. 1] 
expanded to an arbitrary number of customer demand 
distributions. 
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D. SUMMARY AND CONCLUSIONS 



This chapter has shown that the monovarible approxima- 
tion formula provides a numerically exact value for the 
expected shortage cost when the aggregate demand process is 
the sum of arbitrary number of mutually independent Poisson 
processes. This formula utilizes the closure theorem for 
the Poisson distribution, which says that tne sum of inde- 
pendent Poisson random variables is again a Poisson random 
variable. 

Feller [Ref. 5 p.142] shows that the binomial distribu- 
tion with parameters (n,p) converges to the Poisson distri- 
bution with parameter / a= np as n-* co and p 0 in such a way 
that M- =np remains bounded. So the Poisson distribution 
corresponds to a binomial distribution with large n and 
small p. 

Since the Poisson is a limit of binomials and since the 
formula for expected stockout cost is exact for the Poisson 
case, it seems reasonable that the formula might also be 
valid for the case in which demands from each customer 
during a period are binomially distributed. Unfortunately, 
a single counterexample shows that the shortage cost expres- 
sion is not exact for every case in which the demands are 
binomial. To see this, consider the case in which there are 
two customers as follows: 

#Repair cycles Probability of stockout 

customer per period demand per cycle cost 

11 0.9 $ 100 

22 0.1 $ 1,000 

We will show that expression (2. 1) is not valid for this 
example. 

Let X;. be the total demand in the period by customer i . 
Suppose X=X i +X 1 = 2. Then 
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P(X 1 + X i =2) 



k = 0 



?(X x +X 1 = 2 | X x =k) P ( X j_ = K ) 



= 'ZZ ?(X i =2-k)P(X 1 =k) 



= P (X =2) P ( X =0) + P(X =1)P(X =1) 

2 1 2 1 



= P. 



(1“P ) + (t) P (1-p ) 



2 1 



2 1 
= 0.1630 

Now, the conditional probabilities are given by: 



P(X ± = 0, X a = 2) 
P (X =0 | X=2) = — i-1— i — 
i P (X=2) 



.L ( 2*:?Lfi2i:£L = 

P (X = 2 ) 0. 1630 



= 0. 0061 



n#v , ^ P(X 1 = 1)P(X 1 = 1) ftJPin-PilP* nnon 

P (X =1 X=2) = = = 0.9939 

i P (X=2) 0. 1630 

and, the conditional expectation is: 



E[X 1 |X=2] = 0P(X 1 =0|X=2) + IP (X x = 1 1 X =2) = 0.9939 (2.2) 



Now let us compare (2.2) to the value 
sion (2.1); /i- = n/p t ', and E[X 1 |X = 2] = 

1.6364. This is significantly larger than 

( 2 . 2 ) . 



given by expres- 

2 ( = 

V ±(c.d) i - ZCo 1)1 

the exact value 



Thus we can see from the simple example that, for small 
n and appropriate choices of the p^’s, expression (2.1) 
provides a poor approximation for E[ X<; | X=x ] when the total 
demand in a period from each customer is binomially distrib- 
uted. Conseguently , the monovariable expression for the 
total expected shortage cost will also be in error for tnis 
case. However, this example is extreme, and it might be the 
case that for some sets of values for the parameters the 
monovariable expression provides a reasonable approximation 
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for the true shortage cost. We evaluate this possibility in 
a later chapter by simulating actual iemands and shortages 
and comparing the monovaciable expression to the average 
stockout costs computed over a large number of replications. 

Before we proceed to the development of alternative 
approximation formulas and to the simulation results, let us 
first return to the simple example to understand another 
problem that results from the binomial demand situation. 
Suppose there is a single spare part to support both 
customers; i.e. Y=1. As above, suppose that X=2. Mow, what 
is the shortage cost? This cost is either 3 130 or $ 1,300 
depending on the seguence in which the demands occurred 
during the quarter. Below are three possibilities for which 
X=2 ( S corresponds to a demand during the period and ? 
corresponds to no demand) . 



S 

Case 1: 1 | | 

S F 

2 j | | 



p 

Case 3: 2 | 

S S 



S 

Case 2: 1 | j 



In case 3, the stockout cost is clearly $ 1,033. However, 
in cases 1 and 2 the stockout cost could be either $ 100 or 
$ 1,000 depending on when during a repair cycle that the 
demands occurred. If, for example, it is assumed that 
demands occur at the end of a repair cycle, the stockout 
cost in case 1 is $ 100, but in case 2 it would be either 
$ 100 or $ 1,000. If demands are assumed to occur at the 
beginning of a repair cycle, the stockout cost in case 2 
would be $ 1,000, but in case 1 it would be either $ 100 or 
$ 1,000 depending on what kind of rule is used to determine 
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III. derivation of other approximation for solas 

A. INTRODUCTION 

Before simulations are introduced in the next chapter to 
evaluate the accuracy of the mono v aciable formula for the 
binomial "real world" NARF situation, two further approxima- 
tion formulas will be derived. Their accuracy will be eval- 
uated along with the monovariable formula. 

B. DERIVATION OF THE HYPERGEOMETRIC APPROXIMATION FORMULA 

Given two discrete mutually independent demand distribu- 
tions PJX^) and P^Xi) the probability for each possible 
demand vector <x 1 ,x 1 > can be evaluated from the product of 
their probability mass functions, p^x^) P t (x t ) . 

Next in the conditional setting, i.e. given a realiza- 
tion <Xj_,x^> with (x^ + Xi) > Y an interpretation of random 
drawings without replacements suggests an approximation 
based on the hypergeometric distribution. Let 

x = number of balls in urn (=x x *-x 1 ) , 

x. ± = number of balls with associated unit shortage cost , 
x z = number of balls with associated unit shortage cost , 
s = number of shortages [= [x x +x L ) - y} . 

If s balls are drawn out of an urn containing x balls, 
Xj_ of which have associated unit shortage cost n x and x^ of 
which have associated unit shortage cost 7T Z then, according 
to the formula for the expected value of the hypergeometric 
distribution, the expected number of x^^ balls in the sample 
is 



x + x , 
1 2 
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and the expected number of balls is 
x -> 



x + x 
1 2 



shortage cost 



(3.1) 



Multiplying each term by the associated unit 
and summing them yields: 

Xi _ Xi 



S 7Z Z + s 

X + X X + X 

12 12 



This expression represents an approximation for the expected 
shortage costs per period given x ± , x^, Y and the constraint 
that x 1 +x 1 >Y. 

Because and Xj. are mutually independent random vari- 
ables, the unconditional approximation for the expected 
shortage costs per period can be computed by multiplying 
each conditioned approximation with its associated prob- 
ability of occurrence and summing these products, i.e.: 

ri,_ 

J 7”{^*(o/.vxj-y 

y&OXfO 

This sum gives an approximation for the expected shortage 
cost per period. Unfortunately, it does not reduce to a 
simple formula. The hypergeometric formula is programmed as 
the subroutine HYPGEO in the program provided as Appendix A. 
This formula corresponds to the first step of McMasters' 
development of his approximation formula [Ref. 1]. 



Xj^+Xj 






C. DERIVATION OF THE CDF APPROXIMATION FORMULA 

In the conditional setting the hyper geometric approxima- 
tion formula does not reflect any information provided by 
the individual demand distributions. To incorporate this 
information, an additional weight can be introduced into the 
averaging operation in equation (3.1). 
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If (x 1 +x> ) > Y , the quantity Y is split into two 

subquantities Y^ and Y>_ such that x^ /x, = Y^/Y^ . Suppose 

next that the density functions of X ± and X z look like 
Figure 1. The areas represent the average height of the 
density functions in the indicated regions. 

In Figure 1 a relatively small fraction of the prob- 
ability mass over the region [0 / x i ] is located in the region 
(Y J _,x J _] whereas a relatively, large fraction of the prob- 
ability mass over the region [0 / x 1 ] is located in the region 

(^ 2 . > ^ 4 . 3 * 

The motivation for the CDF approximation formula is that 
the low fraction of probability over (Y 1 ,x i ] would result in 
fewer shortages for customer 1 and the high fraction over 
would result in more shortages for customer 2. 
Those fractions of probabilities over the regions (Y ± ,Xj_] 
and (Yi/Xi.] ace incorporated as additional weights into the 
hypergeometric approximation formula. 
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Let CDF1 and C DF2 denote the cumulative distribution 
functions for X x and X x respectively. Then the additional 
weights f ± and £j_ can be computed by 

CDF 1 (x x ) - CDF1 (Y x ) 

1 CDF 1 (xj 

CDF2 (x z ) - CDF2 (Y z ) 

2 CDF2 {x ) 

These weights are then introduced as multipliers in the 
averaging operation of the hyper geometric approximation 
formula, i.e. equation (3.1) becomes 



s 



+ x 



7T ( 



f + 

l 




(3.3) 



12 12 
The expansion to the unconditioned setting is identical 
to the procedure applied when deriving the hypergeometric 
formula at the end of the preceding section (i.e. Equation 
(3.2) ). The CDF approximation formula is programmed as the 
subroutine SOBCDF in the program provided as Appendix A. 



D. COMMENTS QN THE TWO DERIVED APPROXIMATION FORMULAS 

In case of the monovariable formula there is only one 
convex combination of the unit shortage costs computed, 
whereas by applying the hypergeonetric or CDF formulas a 
convex combination of the unit shortage costs for each 
possible demand vector must be computed. It might be that 
computing the convex combinations within the conditioned 
setting provides better approximations than the monovariable 
formula. 

In comparing these approximations under the Poisson 
distribution of Chapter II, the monovariable formula and the 
hypergeometric formula provide the same numerical results in 
the cases of mutually independent Poisson and binomially 
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distributed demands when the latter is such that p, =p for 
all i . 

The suggested CDF formula attempts to refine the hyper- 
geometric approximation formula. However, it requires the 
computation of the sutguantities Yj_ and Yi_ . If Y x and Y z 
are not integers, a interpolation of the individual cumula- 
tive distribution functions has to be performed in order to 
compute the weights and f^. An alternative approach would 
be to round Y ± and Y^ to the nearest integer. Unfortunately, 
both approaches cause distortions, therefore the CDF formula 
was not evaluated by running simulations in the Poisson 
situation. However, it is intuitively appealing and should 
be evaluated along with the monovariable and hypergeometric 
approximation formulas. The main question to be answered as 
far as the CDF approximation formula is concerned, is 
whether the incorporation of the distributions into the 
averaging operation is a step in the right direction. 

E. SDHHARI 

This chapter provided the arguments leading to the two 
further approximation formulas, called the hypergeometric 
approximation formula and the CDF approximation formula. 
Simulation is required to determine the accuracy to be 
expected from applying the three approximation formulas in 
the binomial "real world" situation of the NARF model. 
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IV. DESCRIPTION AND ANALYSIS OF SIMULATIONS 



A. PRELIMINARY REMARKS 

If the individual customer demands are binomially 
distributed with possibly different P[*s and ni’s there is, 
based upon the knowledge of the author, no closed form 
expression available which provides the exact value for the 
expected shortage costs per period. Therefore simulation 
was used in order to estimate the accuracy of the approxima- 
tion formulas . 2 However, as with any simulation, it will 
provide only an estimate of the actual value of the expected 
shortage costs. The precision of the estimate gained from 
simulation is almost exclusively determined by the number of 
simulated periods (replications) . 

B. PROGRAMMING LANGUAGE AND RANDOM NUMBER GENERATOR 

The simulation programs were written in FORTRAN WATFIV. 
WATFIV allows straightforward structured programming which 
improves the readability and verifiability of a program. 
The WATFIV compiler provides a large set of error checks and 
gives, in case of errors, uncoded and precise error 
messages. A further reason is the availability of a WATFIV 
specific debugging mode at the Naval Postgraduate School 
which eases testing of a program. 



2 In the early stages of the research done for this 
thesis the author had written simulation programs for up to 
four customers whose demands per period were assumed to be 
Poisson distributed. The simulations suggested that the 
values of the approximation formula were in fact the exact 
expected shortage costs per period. Driven by this evidence 
the author finally succeeded in finding an analytical proof 
for the case of a Poisson setting. 
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The basic prerequisite for any Monte Carlo type of simu- 
lation is the availability of proper random numbers. The 
simulation programs use 1 LLR ANDOMII* , the Naval Postgraduate 
School Random Number Generator Package [Ref. 6], This 
package has been tested thoroughly and is well validated. 
It generates very good pseudo random numbers. 

C. SIMULATION OUTPUT ANALYSIS FOR MEAN ESTIMATE 

During simulation many in de p ende nt observations (repli- 
cations) of the shortage cost per period are obtained. 
Because, in the context of this thesis, the interest is 
limited to the mean of the shortage costs per period, a 
straightforward output analysis can be applied. 3 

The sample mean is used to estimate the mean of the 
shortage cost per period. It is unbiased and has minimum 
variance. Dividing the sample variance by the number of 
replications gives an unbiased estimate for the variance of 
the mean estimate. It represents the estimated precision of 
the mean estimate derived from simulation. 

1 . Notati on 

C - scalar random variable, shortage cost per period 

C - sample average 

n - number of replications 

sd - standard deviation 

Var- variance 

A 

A hat indicates an estimator, i. e. Var (C) is used 
to estimate Var(C). 

n 

C.-l- Zlc. 

n i= 1 



3 The reasoning presented in this section is mainlv 
extracted from [Ref. 7] as far as output analysis is 
concerned. 
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1 



n 



(C; - C) 2 



Var (C) 



n 



1 



■'T' 

/ 

i=1 



Var (C) 



sd (C) 





Var (C) 

V 

A 

Var (C) 




By the Central Limit Theorem the distribution of the 
sample average of any random variable with a well-defined 
mean and variance approaches more and more the normal 
distribution as the size of the sample becomes larger and 
larger. In general if there are more than twenty indepen- 
dent observations in a sample the normal distribution is 
assumed to provide a sufficient approximation. Therefore, 
when n is sufficiently large, the random variables 

C - E[C] (n - 1) Var (C) 

sd (C) , Var (C) 



are distributed as a standard normal and a chi-sguared 
distribution with (n-1) degrees of freedom, respectively. 

Both statistics are independent, therefore 



"sd (C) 

has a Student's t-distribution with (n-1) degrees of free- 
dom. Because the t-distribution converges to the normal 
distribution with increasing n it can be assumed that 



C - E[C] 
sd (C) 

is distributed as a standard normal random variable. 
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For the unknown true mean, a confidence interval can be 
computed with (1-<rO as the preassigned confidence level: 

P { E[C] within C 1 (^Jl - 1 ' cK 

x. 

where Z 2 represents the standard normal deviate. 

i. 

To increase the precision of the estimated mean of the 
shortage costs per period by a factor of two (i.e. reduce 
the estimated standard deviation of the mean estimate to 
half of its value) the number of replications has to be 
increased fourfold. In general, as far as mean estimation 
through simulation is concerned, the precision of the esti- 
mate is directly proportional to the square root of the 
factor by which the number of replications gets increased. 
All simulations run in connection with this thesis generated 
sample sizes far beyond one thousand. 

D. OUTPUT ANALYSIS FOB SIMULATED TWO BINOMIAL DEMAND 

SOURCES 

There are four different categories of • parameters: 
number of repair cycles per period, probability of one 
demand during one repair cycle, unit shortage cost, and the 
purchased quantity. Each of these parameters is a candidate 
for causing deviations of the numerical results of the three 
approximation formulas from the actual shortage costs per 
period. The tests of all three approximation formulas were 
restricted to two binomial demand sources in order to more 
easily identify those parameters which most affect the accu- 
racy of the approximations. 

The output tableau of the simulation program and expla- 
nations of its entries are provided as Appendix F. The 
simulation results presented in Chapter V for the different 
settings display only those values from this tableau which 
are of intermediate interest. 



28 



V. SIMULATION RESULTS FOR TWO DEMAND SOURCES 



A. PRELIMINARY REMARKS 

The program to perform the simulations for two demand 
sources is provided as Appendix A to this thesis. The 
following sections of this chapter have the following 
str uctur e: 

1. Purpose of the chosen parameter setting; 

2. Presentation of the simulation results; 

3. Analysis of the simulation results. 

B. FIRST SIMULATION SETTING 

In the first simulation setting the number of repair 
cycles and the initial seed are fixed for all runs. The p’s 
and Y vary to analyze their influence on the precision of 
the approximation formulas. In Tables 1-4 the unit shortage 
costs are $ 100 and $ 1,000 for demand source 1 and demand 
source 2, respectively. In Table 5 the unit shortage cost 
for source 2 was reduced from $ 1,000 to $ 200 in order to 
analyze the effect of a lower variability among the unit 
shortage costs on the accuracy of the approximations. 

For each table the purchased quantity Y is a specific 
multiple of the sum of expected demands per period from 
both demand sources. The value of Y is always rounded to 
the nearest integer. However, Y was increased by at least 
one unit in the sequence of Tables 1-4 in order to examine 
the influence of a change in Y on the accuracy of the 
approximations. 

The simulation results presented in this and the 
following section are arranged in tableaus consisting of 
blocks of six numbers which refer to specific parameter 
values. 
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shortage- 4 

cost SD/MEAN s y 

DE VMON j DE VHYP 6 D E VC D F_± 

The meanings and units of the abbreviations are the same 
as described in Chapter IV, section D. In the tableaus. Pi 
and P2 denote the p values for the binomial distributions 
for demand source 1 and demand source 2, respectively. The 
maximal sum of both demands per period is 9 units for 
Tables 1-5. 



1 - Param et er Setting f or T ab les 1- 4: 7 

Number of time units 2 0 

Number of repair cycles (1) 4 

Number of repair cycles(2) 5 

Unit shortage cost(1) $ 100.- 

Unit shortage cost (2) $ 1000.- 

Initial seed for all runs: 8 150449 



4 The observed sample mean for shortage costs per 
period derived from the simulation. 

5 The observed standard deviation for the sample mean 
divided by the sample mean provides a measure (in percent) 
of the precision of the simulation. 

6 Deviations of the values of the three different 
approximation formulas from the simulation result in 
percent . 

7 In order to avoid having very large values for Y in 
Tables 2-4 for the case of P 1 = 0 . 1 and P2 = 0.1, Y = 0 was chosen 
for that case in Table 1. 

8 For all simulations in this section the same initial 
seed was used in order to reduce the variability of the 
simulation results. Differences in the performance of the 
approximation formulas between the tables are therefore 
mostly due to different levels of Y. 



30 



2. Table 1: Y About 157- Below Mean Demand 
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3 . Tab le 2: Y About Equal to Mea n Demand 
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4 . Table 3: Y A bout 15% Above Mean Demand 9 
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9 A string of three asterisks (= 1 ***') indicates , that 
because Y is equal to the maximal demand per period no 
shortages were observed for that specific parameter set. 
Therefore no deviations could be computed. 
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5 . Table 4 : Y A bou t 30% Ab ov e Mean Dema nd 
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6 . Tab le 5: Al l P arameters but Unit Shortage Co s t(2) 
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100,000 replications 



l °The number of replications had to be increased from 
10,000 to 50,000 and finally to 100,000 in order to get the 
same precision level for the simulation result for Tables 

1 - 5 . 
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7 . Analy sis of the Simulation Results: 

The following summary contains the average absolute 
deviation (in percent) of the approximation formulas from 
the simulation results: 







level 
Of Y 


monovar . 
formula 


hyperqeom . 
formula 


CDF 

formula 


Table 


1: 


0.85 mean 


7. 4 


5. 3 


13. 8 


Table 


2: 


mean 


15. 6 


11.4 


19. 5 


Table 


3; 


1.15 mean 


23. 2 


18.5 


24. 5 


Table 


4: 


1.30 mean 


30. 5 


24. 7 


38. 5 


Table 


5: 


1.30 mean 


14. 8 


12.2 


19.0 



with reduced unit shor- 
tage cost variability 

In this setting the hypergeometric formula shows 
slightly more accurate results than the monovariable 
formula. The CDF formula shows significantly less accurate 
results than the two other formulas. 

Tables 1-4 reveal that for all three approximations the 
accuracy of the approximations decreases as I increases. 

The accuracies of the approximations are highly depen- 
dent on the ratio of the two unit shortage costs of the two 
demand sources. If both unit shortage costs were the same, 
the problem is reduced to approximating the expected number 
of shortages per period. The higher the variability of the 
unit shortage costs, the less accurate is the approximation 
of the expected shortage costs per period. To illustrate 
this property, the variability of the unit shortage costs 
was set lower in Table 5 than in Tables 1-4. As can be seen 
whan Table 5 is compared to Table 4, the change of unit 
shortage cost (2) from $ 1,000 to $ 200 reduces the inaccura- 
cies of the approximations to half of their former size. 
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Although Table 1 represents simulation results with a 
high variability of the two unit snortage costs, the average 
absolute deviation for the sixteen runs is 7.4 percent with 
a maximal absolute deviation of 16 percent for the monovari- 
able approximation. For the hyper geometric formula the 

average absolute deviation is 5.3 percent with a maximum 
absolute deviation of 13 percent. These results suggest 
that the monovariable and the hypergeometric formulas 
provide an acceptable accuracy when Y is less than the mean 
of the sum of both demands. Therefore, further analyses 
will be restricted to settings where Y is egual or greater 
than the mean of the sum of demands. 

C. SECOND SIMULATION SETTING 

In this setting the number of repair cycles per demand 
source were egual. To analyze the effect of smaller 
increases of Y the three levels are now: Y= mean, Y= 1.10 

mean (i.e. only 10 percent increase), and Y= 1.25 mean. The 
values of PI and P2 were taken to be 0.1, 0.3, 0.5 and 0.9. 

1 * Par ameter Se tti ng for Tables 6-8: 

Number of time units 10 

Number of repair cycles (1) 10 

Number of repair cycles (2) 10 

Unit shortage cost(1) $ 1000.- 

Unit shortage cost(2) $ 100.- 

Initial seed for all runs: 40150 
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2 . Table_6 : Y E qual to Mean Demani 
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3. Ta ble 7: Y About 10 7 Above Mean Demand 11 
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4 . Table 8: Y A bout 25 % Ab ov e Mean Demand 
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iiThe Y's of Table 7 are at least one unit larger than 
in Table 6: the Y's of Table 8 are at least one unit larger 
than in Table 7. 
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5 . Ana lysis of the Simulation Results: 

The average absolute deviations (in percent) of the 
approximation formulas from the simulation results are: 







level 
of y 


monovar . 
formula 


hyper geon. 
formula 


CDF 

formula 


Table 


6: 


Y = mean 


4. 9 


3.6 


16. 0 


Table 


7: 


Y = 1.10 mean 


7. 5 


3. 4 


17. 1 


Table 


8: 


Y = 1.25 mean 


11.1 


4.4 


26. 2 



Tables 6-8 indicate again the relatively strong depen- 
dence of the accuracy of the approximations on the size of Y 
relative to the mean. Although the variability of the unit 
shortage costs is high, the monovariable and the hypergeome- 
tric formulas seem to provide acceptaole accuracy even for 
the case where Y =1.25 mean. Within this very special 
setting all three approximations indicate tnat the accuracy 
decreases slightly as the difference between Pi and P2 
increases. However, Tables 1-8 do not indicate in all cases 
that the precision of the approximations increase as the p 
values decrease, as might be expected from the fact that the 
Poisson distribution is the limiting distribution for the 
binomial distribution when n becomes large and p becomes 
small. 

D. THIBD SIMULATION SETTING 

This setting is a variation on the first setting. In 
order to get a more general survey about the performances of 
the approximation formulas the initial seeds and all parame- 
ters but the unit shortage costs are different for each run. 
The unit shortage costs are restricted to $ 100 and $ 1,000. 
One row refers to one simulation run. 
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1 . Table 9: Run s w ith W id ely Varyin g Paramstecs 
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4 


0.2 


0. 7 


1 00 


1000 


3. 4 


4 


999 


0 


-12 


-1 


5 


3 


0.6 


0.4 


100 


1000 


4. 2 


5 


255 


14 


27 


22 


9 


4 


0.2 


0.8 


100 


1 0 00 


5. 0 


8 


100 


236 


136 


37 7 


3 


7 


0.6 


0.1 


1000 


100 


2. 5 


3 


777 


32 


1 0 


35 


6 


5 


•0. 4 


0.9 


100 


1000 


6. 9 


7 


654 


30 


1 4 


59 


7 


2 


0. 8 


0.1 


1000 


100 


5.8 


7 


379 


-2 


-12 


-8 


3 


4 


0.6 


0.4 


1000 


100 


3.4 


5 


111 


84 


52 


1 1 1 


6 


1 


0.7 


0.5 


100 


1000 


4. 7 


5 


7 


98 


118 


15 1 


2 


7 


0.8 


0.6 


1000 


100 


5.8 


6 


825 


245 


216 


355 


3 


2 


0.5 


0. 7 


100 


1000 


2. 9 


3 


900 


32 


20 


33 


7 


4 


0.2 


0.9 


1000 


100 


5.0 


5 


693 


-24 


3 


-50 


5 


9 


0.6 


0.3 


1000 


100 


5.7 


7 


387 


17 


4 


33 


3 


1 1 


0. 1 


0.7 


100 


1 000 


8. 0 


9 


951 


-1 


-4 


-1 


6 


7 


0.4 


0.6 


1000 


100 


6. 6 


8 


291 


22 


34 


15 



Each run: 20,000 replications 

2 . Analysis of the Si mulat io n R esults: 

The average absolute deviations (in percent) are: 
DEV MON = 74.7, DEVHYP = 60.1, and DEVCDF = 117.5. These 
error deviations are relatively large. However, a reduction 
in the variability of the unit shortage cost to $ 100 and 
$ 200 would reduce the average absolute deviations by about 
one-half demonstrated in Tables 4 and 5. 



1 2 Abbr eviations in Table 9: 

N1 - number of repair cycles (1) 

N2 - number of repair cycles (2) 

Pi - p parameter (1) 

P2 - p parameter (2) 

Pil - unit shortage costfl) 

Pi2 - unit shortage cost (2) 

MU - expected sum of demand per period 
Y - purchased quantity 

seed - initial seed used for result in current row 
DEVMOH, DEVHYP, and DEVCDF represent the deviations of 
the approximations from the simulation result in percent. 
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E. FOURTH SIMULATION SETTING 

This simulation setting is similar to both the second 
and third setting. As in the second setting, the number of 
repair cycles is the same for both demand sources. As in 
the third setting, the parameters have row by row comparable 
values. A comparison between corresponding rows in Table 9 
and Table 10 should indicate whether the performances of the 
approximation formulas improve if the lower number of repair 
cycles is increased such that both demand sources have the 
same number of repair cycles. The ratio Y/MU is about the 
same as in Table 9. 



1 . Ta ble 10: Id e ntic al Parameters as Tab le 9; Bu t N1=N 2 



N 1 


N 2 


PI 


P2 


Pil 


Pi2 


MU 


Y 


SEED 


DEVMON 


DEVHYP 


DEVCDF 13 


3 


3 


0.8 


0.6 


1000 


100 


4. 2 


5 


175 


12 


1 


34 


4 


4 


(N 

• 

O 


0.7 


100 


1000 


3. 6 


4 


231 


11 


-3 


14 


5 


5 


0.6 


0.4 


100 


1000 


5.0 


6 


725 


-5 


5 


- 15 


9 


9 


0.2 


CO 

• 

o 


100 


1000 


9 


12 


798 


13 


-7 


28 


7 


7 


0.6 


0 . 1 


1000 


100 


4. 9 


6 


183 


7 


- 1 


4 


6 


6 


0.4 


0. 9 


100 


1 000 


7.8 


8 


436 


. 10 


-3 


31 


7 


7 


0.8 


0 . 1 


1000 


100 


6.3 


7 


123 


6 


- 5 


6 


4 


4 


0.6 


0.4 


1000 


100 


4.0 


5 


123 


9 


0 


17 


6 


6 


0. 7 


0.5 


100 


1 000 


7.2 


8 


190 


-4 


4 


-20 


7 


7 


0.8 


0.6 


1000 


100 


9.8 


1 0 


850 


3 


-2 


21 


3 


3 


0.5 


0. 7 


100 


1000 


3. 6 


4 


136 


7 


- 2 


17 


7 


7 


0.2 


0.9 


10 00 


100 


7. 7 


8 


212 


-17 


15 


-4 1 


9 


9 


0.6 


0.3 


1000 


100 


8. 1 


10 


187 


1 


-7 


7 


11 


1 1 


0 . 1 


0. 7 


100 


1 000 


8. 8 


10 


449 


7 


- 1 


5 


7 


7 


0.4 


0.6 


1000 


100 


o 

• 


8 


725 


-1 


7 


- 1 1 



Each run: 20,000 replications 



13 The abbreviations are the same as for Table 9. 
(see preceding footnote) 
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2 . Analysis of the Simu la tion R e su lts: 

The average absolute deviations (in percent) are: 
DEV HON = 7.5, DEVHYP = 4.2, and DEVCDF = 18.1. The simu- 
lation results establish that the monovariable and the 
hyper geometric approximation formulas provide better accu- 
racy when the number of repair cycles is the same for both 
demand sources. 

F. FIFTH SIMULATION SETTING 

The increase of the lower number of repair cycles of 
Table 9 such that both demand sources had an equal number of 
repair cycles resulted in a dramatic improvement in the 
accuracy of all three approximation formulas. This improve- 
ment could be caused by two facts. The first is that the 
approximations improve when there is an identical number of 
repair cycles, as established in the second and fourth simu- 
lation settings. The second is that an increase of the 
parameter n is one condition for the convergence of the 
binomial distribution to the Poisson distribution. 
Therefore the improvement could be caused simply by the 
increased number of repair cycles. 

Table 11 displays results of simulation runs with the 
number of repair cycles of one demand source of Table 10 
further increased. The ratio Y/MU is about the same as in 
Table 10. Comparison of Tables 9,10, and. 11 provides an 
indication of both the impact of increasing the n values and 
having them equal in value. 
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1 . T ab le 11: Co ntin uation of Table 10; But N1 or N2 

Incre ased 



N 1 


N 2 


PI 


P2 


Pil 


Pi2 


MU 


Y 


SEED 


DEVMON 


DEVHYP 


DEVCD 


if 


3 


5 


0.8 


0.6 


1000 


100 


5.4 


6 


212 


43 


w 


31 


w 


77 


W 


5 


4 


0.2 


0. 7 


1 00 


1000 


3.6 


5 


542 


-2 


I 


-13 


w 


-3 


I 


5 


7 


0.6 


0.4 


100 


1 000 


5. 8 


7 


239 


-13 


w 


-7 


w 


-23 


W 


9 


10 


0.2 


0.8 


100 


1000 


9.8 


13 


822 


-1 1 


I 


-26 


.K 


-2 


I 


9 


7 


0.6 


0. 1 


1000 


100 


6. 1 


7 


17 


3 


I 


-3 


K 


2 


I 


6 


7 


0.4 


0.9 


100 


1000 


8.7 


9 


489 


-7 


I 


-17 


W 


9 


I 


8 


7 


0.8 


0. 1 


1000 


100 


7. 1 


8 


899 


0 


I 


-10 


w 


-2 


I 


4 


6 


0.6 


0.4 


1000 


100 


4.8 


6 


117 


26 


w 


15 


w 


41 


w 


6 


8 


0.7 


0.5 


1 00 


1000 


8.2 


9 


299 


-15 


V 


-10 


w 


- 26 


» 


7 


9 


0.8 


0.6 


1000 


100 


11 


1 1 


168 


17 


w 


1 1 


w 


33 




5 


3 


0.5 


0. 7 


100 


1 000 


4.6 


5 


256 


20 


w 


11 


w 


32 


w 


7 


8 


0.2 


0.9 


1000 


100 


8.6 


9 


437 


8 


I 


51 


w 


-21 


I 


11 


9 


0.6 


0.3 


1000 


100 


9.3 


1 1 


219 


-1 


E 


-7 




1 


X 


22 


1 1 


0. 1 


0.7 


100 


1 000 


9.9 


1 1 


612 


15 


w 


5 


w 


22 


» 


9 


7 


0.4 


0.6 


1000 


100 


7.8 


9 


439 


-13 


w 


-7 


E 


-23 


r.7 

* 



Each run: 20,000 replications 

2 . Anal ysis of the Simulat io n Results: 

The average absolute deviations (in percent) are: 
DEVMON = 12.9, DEVHYP = 14.9, and DEVCDF = 21.1. 

In general, the accuracy of the approximations 

decreased. However, when compared .to Table 10, the monova- 
riable formula improved its performance in 6 out of 15 
cases, and the CDF formula in 7 out of 15 cases, whereas the 
hyper geometric formula lost accuracy in 14 out of 15 cases. 



1 4 The additional letters in columns DEVMON, DEVHYP, and 
DEVCDF indicate: 

I - improved 
E - same level 
W - became worse 
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It can be concluded that an increase of the number repair 
cycles has a larger improving influence on the accuracy of 
the monovariable ana the CDF formulas than it has on the 
accuracy of the hypergeometric formula. The accuracy of the 
hypergeometric formula is more sensitive to the ratio of N1 
to N2 than the two others; if performs best when N1 = H2. 

The impact of increasing n can guickly oe seen when the 
average absolute deviations of Table 11 are compared to 
Table 9. 



DEV 110N DEV HYP 

Table 9: 74.7 60.1 

Table 11; 12.9 14.9 



DEVCDF 
117.5 
21 . 1 



Table 9 and 11 display results for very similar param- 
eter values. In corresponding rows all parameters values are 
identical except N 1 , N2 and Y. However, in each row of 
Table 11 there are increased but unequal numbers of repair 
cycles, the Y's of Table 11 are larger than in Table 9 in 
order to keep the ratio Y/MU at about the same level for 
corresponding rows. 



G. SIXTH SIMULATION SETTING 

Up to this point the timing of demands has been assumed 
to be at the start of the repair cycle. However, this may 
not always be the case at the NARF; the timing is dependent 
on the length of time for teardown and reassembly of a 
repairable. It is quite possible that demand for an item 
could occur at the end of a cycle when the reapairable is in 
its final stage of reassembly. In addition, in Chapter II 
we illustrated the impact of the timing of the demands 
during the repair cycle on the shortage costs with a simple 
example. 
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Another assumption has been that, in the event of two 
customers simultaneously demanding the last unit in inven- 
tory , the supply department would toss a coin to see who 
gets the item. However, it is not unreasonable that the 
supply department might favor the customer with the higher 
expected demand per period. 

In this simulation setting demands are assumed to occur 
at the end of the repair cycle and, in the case of a tie for 
the last unit, the probability that a customer will get it 
is the ratio of his mean demand for the period to the total 
mean demand from all customers. The simulations were reruns 
of Table 9 with "Stage 1" involving the timing change and 
"Stage 2" involving both the timing change and the new 
priority rule. 

1 . Ana ly sis of the S imulat io n Results: 

The average absolute deviations (in percent) are: 









DEV HON 


DEV HYP 




DEVCDF 


In 


Table 


9: 


74.7 


60. 1 




117.5 


For 


Stage 


1 : 


17.8 


18.9 




21.5 


For 


Stage 


2: 


10. 1 


17. 1 


# 


17.7 



These results establish with very high significance the 
decisive influence of the timing of the demands on the accu- 
racies of ail three appproximation formulas. They also 
suggest that, if the supply department would consider a 
priority scheme based on expected demands rather than 
tossing a coin for the last unit, the shortage cost approxi- 
mations could be used with more confidence. 

Interestingly, the occurrence of demand within the 
repair cycle does not matter, when the number of repair 
cycles is identical because the cycles are then of identical 
length too. In the case of an equal number of repair 
cycles, the order in which demands occur will not change as 
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2 . 



Table 12: Re run of Table 9 with Changed "Real W orld " 



Accuracies 
achieved 
in Table 9 



Stage 1, 
shifted 
occurrence 
of demand 



Stage 2, 
additionally 
changed 
priority rule 



'-MON 


D-HYP 


D-CDF | 


D-MON 


D-HYP 


D-CDF | 


D-MON 


D- HYP 


d-cd: 


283 


230 


506 


-31 


-40 


10 


2 


- 13 


6 1 


0 


-12 


-1 


16 


2 


15 


5 


8 


4 


14 


27 


22 


-21 


-12 


-15 


-9 


2 


-3 


236 


136 


377 


25 


- 13 


77 


15 


-19 


64 


32 


10 


36 


4 


-13 


7 


1 


- 15 


5 


30 


14 


59 


10 


-4 


34 


5 


-9 


28 


-2 


-12 


-8 


33 


20 


26 


5 


-6 


- 1 


84 


62 


1 1 1 


-1 


-13 


14 


-4 


-15 


1 1 


98 


118 


151 


-61 


-57 


-50 


-45 


-40 


-31 


245 


216 


355 


-37 


-42 


-17 


-30 


-35 


-7 


32 


20 


38 


-9 


-17 


-5 


-8 


- 16 


-4 


-24 


3 


-50 


-1 


34 


-35 


7 


45 


-30 


17 


4 


33 


-6 


-16 


8 


-6 


- 16 


7 


-1 


-4 


- 1 


11 


8 


10 


7 


4 


7 


22 


34 


15 


-2 


7 


-9 


3 


13 


-3 



Each run: 20,000 replications 



long as demands occur at the same point in time within the 
repair cycles for both demand sources. This fact explains 
the relatively good performance of the approximation 
formulas in simulation settings 2 and 4. In the case of an 
unequal number of repair cycles the order will generally 
change. 

H. SEVENTH SIMULATION SETTING 

As the number of repair cycles increases, the influence 
of the timing of demands within a repair cycle on the accu- 
racy of the approximations will decrease. 
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Simulation settings 1-6 provided strong evidence that 
the CDF approximation formula is clearly the worst. It 
does not provide sufficient accuracy and will be excluded 
from further simulations and analysis. 

The following restrictions for the parameter values 
should ensure an acceptable accuracy for the monovariable 
and the hypergeometric approximation formulas when they are 
applied to the binomial NARF model. The restrictions are 
derived from simulation settings 1-6: 

1. The number of repair cycles for each demand source is 
greater than or egual to 10. 

2. The mean of the sum of both demands per period plus 
1.50 times its standard deviation is an upper bound 
for Y. is 

3. The ratio of the larger unit shortage cost to the 
smaller is at most three. 

Table 13 represents results of simulations which were 
run with parameter values within the above restrictions. 

1 . Analy sis of the Si mu lati on Res ul ts: 

The average absolute deviations (in percent) are: 
DE7M0N = 6.3 and DEVHYP = 5.8. The maximal absolute devia- 
tion is 14 percent for the monovariable formula and 
16 percent for the hypergeometric formula. 

Because the fifteen simulation runs are independent 
observations for DEVMON and DEVH YP , nonparametric statis- 
tical analysis can be performed. Conover [Ref. 8 p.118] 
provides a nonparametric method for computing tolerance 
limits. His method applied to Table 13 establishes that, 
with a confidence level of 95 percent, at least 81 percent 



15 If Y is egual to the mean of the sum of demands per 
period. plus 1.5 of its standard deviations, shortages will 
gccur in about one out of 15 periods. This quantification 
is a rough estimate derived from the standardized normal 
distribution, i. e. P (Z>1 . 5) =0. 067. 
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2 • Table 13: Run s wi th Re s tr icted Varyi ng 
Parameters 1 6 



N 1 


N 2 


PI 


P2 


MU 


STD 


Y 


Pil 


Pi2 


Seed 


SHCOST 


DEVMCil 


D E7 HY ? 


10 


15 


0.2 


0. 8 


14.0 


2.0 


16 


100 


200 


909 


28 


-1 


-5 


12 


18 


0.4 


0.6 


15.6 


2.7 


19 


300 


100 


109 


19 


5 


9 


15 


12 


0.7 


0.6 


17.7 


2.5 


20 


100 


300 


89 1 


33 


14 


16 


10 


20 


0.3 


0. 5 


13.0 


2.7 


17 


300 


100 


510 


11 


6 


1 1 


1 1 


22 


0 . 1 


0. 9 


20.9 


1.7 


23 


100 


300 


175 


21 


3 


-2 


21 


1 4 


0.3 


0.2 


9 . 1 


2.6 


12 


300 


100 


175 


43 


0 


- 1 


12 


16 


0.8 


0 . 1 


1 1.2 


1.8 


13 


100 


250 


311 


20 


-12 


0 


20 


40 


0. 9 


0.9 


54.0 


2.3 


56 


1 00 


300 


1 1 3 


55 


-13 


-1 3 


10 


1 5 


0.2 


0.2 


5.0 


2.0 


8 


100 


300 


975 


16 


-2 


-2 


20 


1 5 


0.2 


0.6 


13.0 


2.6 


16 


300 


100 


305 


30 


-12 


-5 


12 


1 8 


0. 9 


0.5 


19.8 


2.4 


23 


1 00 


300 


129 


18 


- 1 1 


-4 


13 


39 


0.4 


0.7 


32 .5 


3.4 


37 


1 00 


300 


610 


36 


-5 


-7 


11 


33 


.05 


0 . 1 


3.9 


1.9 


7 


300 


100 


7 1 1 


6 


5 


7 


18 


12 


0. 2 


0.3 


7.2 


2.3 


8 


300 


100 


1 1 1 


1 1 8 


-2 


0 


10 


40 


0.3 


0.5 


23.0 


3.5 


25 


100 


300 


637 


173 


-4 


-5 



Each run: 20,000 replications 



of the absolute deviations of the monovariable and hypergeo- 
metric formulas will not exceed the maximal absolute devia- 
tions of Table 1 3. 

I. SUMMARY AND CONCLUSIONS 

This chapter has shown that the monovariable and the 
hypergeometric formulas provide significantly more accurate 
results than the CDF formula does. The following additional 



1 6 STD denotes the standard deviation for the sum of botn 
demands per period. SHCOST denotes the average value for 
the shortage costs per period as observed during simulation. 
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remarks refer to the monovariable and the hy pergeomet ric 
formulas. 

All simulations show a clearly recognizable trend that 
the accuracy decreases as Y increases with levels of Y which 
correspond to most "real world" situations (i.e. mean - 15 %, 
mean +30%) . This suggests that the approximations are not 
good in the upper tails of the demand distributions. 

The simulations suggested that the approximations 
provide their most accurate results when there is an iden- 
tical number of repair cycles. Even in the case of a high 
variability of the unit shortage costs (i.e. 3 1,000 and 
$ 100 ), the average absolute deviation from the simulation 
results ( simulation settings 2 and 4) is about 10 % for the 
monovariable formula and 5 % for the hypergeometric formula. 

The variability of the unit shortage costs is an impor- 
tant factor for the accuracy of the approximation formulas . 
Reducing the variability of the unit shortage costs from 
$ 1,000 and $ 100 to $ 200 and 5 100 almost doubles the 
accuracy to be expected from both approximations, as was 
demonstrated in simulation setting 1. 

Simulation settings 4 and 6 suggest that the approxima- 
tions perform significantly better when the demands are 
assumed to occur at the end of the repair cycle. The influ- 
ence of the timing of the demands on the accuracy of the 
approximation formulas decreases as the number of repair 
cycles increases. 

In settings with different numbers of repair cycles per 
period the monovariable and the hypergeometric formulas 
provide about the same level of accuracy. By imposing the 
following restrictions for most practical applications an 
acceptable level of accuracy can be expected: 

1. The number of repair cycles for each demand source is 
greater than or equal to 10. 



46 



2 . 



The mean of the sum of both demands per period plus 
1-50 times its standard deviation is an upper bound 



for Y. 

3. The ratio of the larger unit shortage cost to the 
smaller is at most three. 

More specifically, when observing these restrictions it 
is to be expected that, with a confidence level of 95 
percent, at least 81 percent of the actual absolute devia- 
tions of the monovariable formula will not exceed 14 
percent, and, in the case of the hypergeometric formula, 
they will not exceed 16 percent. 



VI. MORE THAN TWO BINOMIAL DEMAND SOURCES 

A. PRELIMINARY REMARKS 

This chapter analyzes the situation when there are more 
than two demand sources and presents results from the asso- 
ciated simulation runs. In particular, it considers the 
approximation of the distribution of the sum of demands per 
period by the normal distribution. 

The computational effort required to use the hypergeoae- 
tric formula increases very rapidly with an increasing 
number of demand sources. Therefore only the aonovariable 
formula will be evaluated for the case of more than two 
individual demand sources. 

B. METHODS FOR ESTIMATING THE EXPECTED SHORTAGE COSTS PER 

PERIOD 

Two basic ways for estimating the expected shortage 
costs per period in the case of "many" binomial distributed 
demand sources are analyzed in this section. The first is 
to simulate the complete period for a sufficient number of 
replications to provide the desired precision. If the 
desired precision has been prespecified, the number of 
necessary replications can be determined by the program 
during its execution. This simulating approach gives an 
unbiased estimate for the expected shortage costs per period 
and, additionally, allows for varying the time at which the 
demands are assumed to occur within the repair cycle. 

The second way to estimate the expected shortage costs 
pec period is to apply the monovariable approximation 
formula. However, before applying the monovaciable formula 
the distribution for the aggregate demand has to be 
computed. 
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In the case of "many" binomially distributed demand 
sources with different p parameters there are essentially 
four approaches for computing the distribution of aggregate 
demand: 

1. The first approach is to estimate the distribution of 

the sum of demands per period by §1 ®uld.tion . The 

demand occurrences during the period are simulated 
and added up to get one observation for the sum of 
demands per period. This is repeated for a specified 
number of replications. The observed empirical 
distribution for the sum of demands during simulation 
can be taken as an estimate for the distribution of 
the sum of demands per period. 

2. The second approach is to compute the exact distribu- 
tion of the sum of demands. However, the number of 
computations for computing the exact distribution 
grows as a factorial expression whereas the number of 
computations performed during simulation only grows 
approximately linearly along with the number of 
demand sources. For instance, if there are five 
demand distributions with <6,7,8,9,10> repair cycles 
per period, then a five-fold convolution is required 
which incorporates 7x8x9x10x11 = 55,440 iterations. 

3. There are two analytical approximating methods which 

can reduce the huge number of computations involved 
in finding the distribution of the sum of the demands 
per period. One approach is to approximate the indi- 
vidual binomial distributions by Po is so n 

distr ibutions . 

4. The second analytical method is to app ro x imat e the 
individual distributions by nor ma l d ist ribu tions . 
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C. APPROXIMATION OF THE SOM OF DEMANDS BY THE NORMAL 

DISTRIBUTION 

An advantage of the monovariable approximation formula 
is that it does not require the consideration of all 
possible permutations. It requires oniy the distribution of 
the sum of demands per period. Therefore it can be applied 
when the individual binomial distributions are approximated 
by normal distributions. The sum of normal random variables 
is again normally distributed. Its mean and variance are 
the sum of the individual means and variances [Ref. 9 
p. 405 ]. 

The Poisson distribution could also be used as an 
approximation since the sum of Poisson distributed random 
variables is a Poisson random variable whose mean is the sum 
of the individual means. [Ref. 9 p.405 ]. 

The advantage of the normal distribution over the 
Poisson is that it has two parameters. One parameter deter- 
mines the level ( = mean) and the other one the spread 
(=variance) . The effects of the two parameters are mutually 
independent and permit the modelling of a wider variety of 
demand behavior. Therefore the normal distribution is often 
preferred over the Poisson distribution when one attempts to 
fit a theoretical distribution to an empirical distribution. 

A continuity correction was used to reduce the error of 
approximating a discrete distribution with a normal distri- 
bution. If, for instance, the probability associated with 
the integer I is to be computed, the normal approximation 
with continuity correction uses the area under the normal 
density from (I - 0.5) to (I + 0.5) to approximate the 
probability. 

A further reason for using the normal distribution is a 
consequence of the Central Limit Theorem. From the Central 
Limit Theorem it follows that as the number of demand 
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sources increases, the aggregate distribution of demands per 
period will be increasingly better approximated by the 
normal distribution. 

D. RUNNING THE SIMULATION 

In case of an arbitrary number of demand distributions 
the monovariable approximation formula uses an approximating 
distribution to evaluate the expected number of shortages 
per period. Therefore, deviations between the result of 
simulation and the monovariabie approximation formula when 
aggregate demand is assumed to be normally distributed might 
be caused by both the process of applying the approximation 
formula and this use of the normal distribution. 

The output tableau of the simulation program is basi- 
cally an extension of the output tableau described in 
Chapter 4, section D of this thesis. To enable an efficient 
output analysis, the program provides more key statistics 
than the former program for two binomial demand 
distributions. 

NGFDM D= #0F TIME 0NITS= 30 2 REPLI=5000Q I-SEED= 909 



#0F 


REP- 


PROBA- 


UNIT-SHORT. 


FIT 3 


CYCLES 


BILITY 


MEAN COST 


N *P 


N*Q 


1 : 


10 


1: 0.20 


1: 2.00 1: 100.00 


2.0 


8. 0 


2: 


15 


2: 0.80 


2: 12.00 2: . 200.00 


12.0 


3.0 


TOTAL 


MEAN = 


14.00 


TOTAL S- DE V= 2.00* Y= 16 


MAX-DMD= 


25 



SIMULATION: AVE-COST= 28.365 #0F SHOE* S= 0.1525 6 

SHSIM U= 10. 30* SD— COST= 0.42 7 SD-#0F SH'S = 0.0023 7 

NORMAL-APPROXIM. : AVE-COST= 30.01 #0F SH0R’S= 0.1616 

SUMMARY STATISTICS: 

SD/A V E= 1.47%s PISIMA= 185.97$« PINORA= 135.71$*° 

DEVMON= 5.83%** DSVNSH= 5.97S* 2 DEVPIA= -0.14%* 3 
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1. Number of demand sources; 

2. Number of time units per period; 

3. The last two columns indicate how well the individual 
binomial distribution might be approximated by an 
normal distribution. In general, n*p and n* (1-p) 
should both have at least a value of five to ensure a 
sufficient fit for the individual distribution. 
However, applying the monovariable formula incorpo- 
rates only the distribution of aggregate demands. 
With an increasing number of demand sources it is to 
be expected that according to the Central Limit 
Theorem the approximation of the aggregate demand 
distribution will be by far better than indicated by 
the N* P and N* (1 -P) values for the individual 
distributions. 

4. Standard deviation of the sum of demands; 

5. Estimated expected shortage costs per period; 

6. Average number of shortages per period as they were 
observed during simulation. 

7. The observed standard deviations for the estimates 
give a numerical hint about the precision of the 
estimates. 

8. The observed standard deviation of the sample mean 
divided by the sample mean of the shortage costs per 
period. 

9. Average unit shortage costs derived from the simula- 
tion, i.e. accumulated observed shortage costs 
divided by the observed number of shortages. 

10. Average unit shortage costs weighted by the means of 
the individual demand distributions corresponding to 
the way the monovariable approximation formula uses 
these weights. 

11. Deviation of the result of the monovariabie approxi- 
mation formula from the estimate for the expected 
shortage costs per period derived from simulation. 
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12. The deviation of the average number of shortages used 
by the monovariable approximation formula from the 
average number of shortages as observed during 
simulation. 

13. The deviation of the mean weighted average unit 
shortage cost computed when applying the monovariable 
approximation formula from the observed average unit 
shortage cost. DEVNSH and DEVPIA provide a rough 
means for distinguishing whether the deviation 
between simulation results and monovariable approxi- 
mation is caused by the approximation formula or by 
the approximation of the distribution of the aggre- 
gate demands per period by the normal distribution. 

E. SIMULATION RESULTS 

This section presents conclusions derived from simula- 
tion runs performed in order to evaluate the performance of 
the monovariable approximation formula when the total demand 
distribution is approximated by the normal. For all simula- 
tions presented in this section the range of possible param- 
eter values is within the limits of the three restrictions 
listed at the end of Chapter V. 

Table 14 is a direct continuation of Table 13. The only 
difference is that the distribution of the sum of demands as 
now approximated by a normal distribution. 

The average absolute deviations (in percent) are: 
DEVMON = 11.1, DEVNSH = 10.4, and DEVPIA = 5.8. The ’DEVPIA* 
column of Table 14 is almost identical to the 'DEVMON* 
column of Table 13. The 'DEVNSH* column of Table 14 measures 
the additional inaccuracy incurred by applying the normal 
approximation instead of using the exact distribution for 
the sum of demands as was done in Table 13. 
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1 . Table_14 j__Runs_with_I^entical_?a rame ters_as 
Table 13 * 7 



N 1 


N 2 


PI 


P2 


HD 


STD 


Y 


Pil 


Pi2 


Seed 


DEVMON 


DEVNSH 


devp: 


10 


15 


0.2 


0.8 


14.0 


2.0 


16 


100 


200 


909 


5. 8 


6.0 


- 0 . 1 


12 


18 


0.4 


0. 6 


15.6 


2.7 


19 


300 


100 


1 09 


9.4 


5. 4 


3. 8 


15 


12 


0.7 


0.6 


17.7 


2.5 


20 


100 


300 


891 


22. 9 


10.8 


1 1.0 


10 


20 


0.3- 


0.5 


13.0 


2.7 


17 


300 


100 


510 


5. 1 


-1.7 


6. 9 


11 


22 


0.1 


0.9 


20.9 


1.7 


23 


100 


300 


175 


22. 9 


24. 1 


-1.0 


21 


14 


0.3 


0.2 


9 . 1 


2.6 


12 


300 


100 


175 


-8. 1 


-7.7 


-0. 4 


12 


1 6 


0. 8 


0 . 1 


1 1.2 


1.8 


13 


100 


250 


311 


-8.4 


2.7 


- 10. 3 


20 


40 


0. 9 


0.9 


54.0 


2.3 


56 


100 


300 


1 13 


-2. 1 


11.2 


-12. 0 


10 


15 


0.2 


0.2 


5.0 


2.0 


8 


1 00 


300 


975 


-23. 1 


-20.1 


-3. 8 


20 


15 


0. 2 


0.6 


13.0 


2.6 


16 


300 


100 


305 


-14.2 


-6 . 1 


-3. 6 


12 


1 8 


0.9 


0.5 


19.8 


2.4 


23 


100 


300 


129 


-1.0 


14.4 


-13.5 


13 


39 


0.4 


0.7 


32 .5 


3.4 


37 


1 00 


300 


610 


3. 8 


8.8 


-4. 6 


11 


33 


.05 


0 . 1 


3.9 


1.9 


7 


300 


100 


711 


-33. 4 


-36. 1 


4. 2 


18 


12 


0.2 


0.3 


7.2 


2.3 


8 


300 


100 


1 1 1 


-2.3 


0. 2 


-3. 1 


10 


40 


0.3 


0.5 


23.0 


3.5 


25 


100 


300 


637 


-3. 9 


-1.1 


-2. 8 



Each run: 20,000 replications 

If in one row of Table 14 the actual values of DEVNSH 
and DEVPIA are either both negative or both positive then 
the actual overall deviation (DEVMON) represents an accumula- 
tion of DEVNSH and DEVPIA. Otherwise, when the signs of the 
actual values for DEVNSH and DEVPIA are different then the 
aggregate affect of DEVNSH and DEVPIA on the value of DEVMON 
results in a partial compensation of their individual 
affects. In Table 14 the actual values for DEVNSH and 
DEVPIA have different signs in nine out of fifteen rows. 



1 7 DEVNSH = deviation in number of short's (in percent) ; 
DEVPIA = deviation in average unit shortage cost 
(in percent) . 



54 



The average absolute deviations of DEVNSH and DEVPIA 
indicate that about two thirds of the overall average abso- 
lute deviation is caused by using the normal distribution as 
an approximation for the distribution of the sum of demands 
and the rest from applying the monovariable approximation 
formula. It can be concluded that in the case of two demand 
sources it is worthwhile to compute the exact distribution 
of the sum of demands if the monovariable approximation 
formula is going to be applied. 

Simulation results for three demand sources are 
displayed in Appendix C. The mean of the distribution of 
the sum of demands is equal to 21.6 demands per period and 
its standard deviation is equal to 3.75. These tables show 
that the accuracy of the overall approximation decreases 
with increasing Y. The following results (in percentages) 
were obtained: 





DSVMON 


D EV NS H 


DEVPIA 


20 


-1.3 


0.3 


-1.6 


23 


-3.4 


-1.2 


-2.2 


26 


-7. 4 


-4.7 


-3. 1 



These results indicate again that the accuracy decreases 
as Y increases. 

Appendix D contains results of fifteen simulations for 
three demand sources. The fifteen cases of Table 14 were 
expanded by one additional demand source in order to provide 
a basis for comparing the results of Appendix D with Table 
14. For each tableau of Appendix D the ratio Y/M(J is about 
the same as in corresponding rows of Table 14. 

The average absolute deviations (in percent) are: 
DEVMON = 5.8, DEVNSH = 4.2, DEVPIA = 3.8. For Table 14 they 
were: 11.1, 10.4, and 5.8, respectively. The expansion to 
three demands sources shows already a strong numerical 
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influence from the Central Limit Theorem. The average abso- 
lute deviation for the number of shortages dropped from 
10.4% to 4.2%. In six out of fifteen cases the actual 
values of DEVNSH and DEVPIA have different signs which 
reduces the overall deviation for those cases. 

Conover's method [Ref. 8 p.118] for computing tolerance 
limits applied to Appendix D establishes that, with a confi- 
dence level of 95 percent, at least 31 percent of the abso- 
lute deviations (DEVMON) will not exceed 11.9 percent if the 
three restrictions for the ranges of parameter values are 
applied. This accuracy should be acceptable for most 
applications. 

Appendix E contains results of simulation runs for five 
demand sources with one parameter set and three levels of Y, 
i.e. mean, mean plus one standard deviation, and mean plus 
1.5 standard deviations. The deviations (in percent) are: 





DEVMON 


DEVNSH 


DEVPIA 


42 


CT> 

• 

O 


o 

• 

o 


0.9 


46 


0.5 


-1.7 


1.2 


48 


CM 

t 

O 

1 


-2.3 


2.2 



These results indicate that an acceptable level of accu- 
racy can be expected in the cases of five and more demand 
sources . 



F. SUMMARY AND CONCLUSIONS 

This chapter has shown that in case of two demand 
sources it is worthwhile to compute the exact distribution 
of the aggregate demands. In the case of three demand 
sources, the benefit gained from computing the exact distri- 
bution for the sum of demands is significantly lower. 
However, if the numbers of repair cycles are small - which 
means a relatively small computational effort is required in 
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computing the exact distribution of the sum of demands - it 
is recommended that the exact distriDution for the sum of 
demands be computed instead of using the normal approxima- 
tion, Otherwise for the cases of three or more demand 
sources the normal approximation should be applied. 

All simulation results presented in this chapter indi- 
cate that the sign of the actual deviations caused from 
approximating the distribution of the aggregate demands 
(DEVNSH) is independent from the sign of the deviation 
resulting from applying the monovariable approximation 
formula (DEVPIA) . If for a given parameter setting the 
signs of DEVNSH and DEVPIA are different then the overall 
deviation (DEVMON) represents a partial compensation of the 
aggregated affects of DEVNSH and DEVPIA; otherwise, the 
overall deviation is about egual to the sum of DEVNSH and 
DEVPIA. Therefore, the impact of the normal approximation 
on the variability of the accuracy will be stronger than its 
impact on the expected level of accuracy. 
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VII. SUMMARY AND SUGGESTIONS FOR FURTHER ANALYSIS 



A. SUMMARY 

Chapter II proved that the monovariable approximation 
formula for the expected shortage costs per psriod provides 
the exact value for the expected shortage costs per period 
in the case of any number of Poisson distributed individual 
demands. Because the monovariable approximation does not 
provide the exact value for the expected shortage costs per 
period in the case of binomially distributed individual 
demands as they are modelled by tna NARF model, two further 
approximation formulas, the hypergeometric formula and the 
CDF formula, were introduced in Chapter III. 

Chapter IV described the simulation methodology which 
was used in Chapter V to evaluate the performances of the 
three approximation formulas in the case of two binomially 
distributed demand sources. Chapter V showed that the mono- 
variable and the hypergeometric formulas provide about the 
same level of accuracy. Both formulas provide significantly 
more accurate results than the CDF formula. They are there- 
fore the ones that should be used in applications. 

Chapter V showed that three factors are significant in 
determining the accuracy of the three approximation 
formulas: 

1. The occurrence of demand at the end of the repair 
cycle improves the accuracy of the approximations. 
However, in almost all "real world" situations which 
are modelled by the NARF model, demand will occur at 
the beginning of the repair cycle. The approximation 
error incurred hereby will decrease as the numbers of 
repair cycles increase. 
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The accuracy of the approximations for the expected 
shortage costs per period decreases as the available 
quantity Y increases reflecting the fact that the 
approximations are not so good in the upper tails of 
the distributions. 

3. If both unit shortage costs were the same, the 
problem of approximating the expected shortage costs 
per period is reduced to approximating the expected 
number of shortages per period. The higher the vari- 
ability of the unit shortage costs, the less accurate 
is the approximation for the expected shortage costs 
per period. 

By imposing three restrictions on the range of parameter 
values, an acceptable expected level of accuracy can be 
achieved for the monovariable and the hy pergeometr ic formula 
in the cases of two demand sources. These restrictions are: 

1. The number of repair cycles for each demand source is 
greater than or equal to 10. 

2. The mean of the sum of both demands per period plus 
1.50 times its standard deviation should be an upper 
bound for Y. 18 

3. The ratio of the larger unit shortage cost to the 
smaller is at most three. 

Tolerance limits for the accuracy were computed from 
Simula tion results in Chapter V. With a confidence level of 
95 percent, at least 81 percent of the actual results from 
the monovariable approximation formula will incur at most an 
absolute deviation of 14 percent from the exact numerical 
value of the expected shortage costs per period, and at most 
16 percent in case of the hypergeometric formula. As a 



18 If Y is equal to the mean of the sum of demands per 
period plus 1.5 of its standard deviations shortages will 
occur in about one out of 1 5 periods. This quantification 
is a rough estimate derived from the standardized normal 
distribution, i. e. P (Z> 1 . 5) =0. 067. 
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general tendency, the hy pergeometric formula seems to 
perform slightly better than the monovariable approximation 
formula, especially in the case when there are equal numbers 
of repair cycles associated with each demand source. 

Chapter VI expanded the setting to more than two binomi- 
ally distributed demand sources. The problem with the 
hypergeometric formula is that all possible pec mutations of 
demands must be evaluated whereas the monovariable formula 
requires only the distribution of the aggregate demands per 
period. Because the number of computations needed to apply 
the hypergeometric formula increases rapidly as the number 
of demand sources is increased, only the monovariable 
formula was evaluated in the cases of more than two demand 
sources . 

Chapter VI showed that in case of two demand sources it 
is worthwhile to compute the exact distribution of the 
aggregate demands rather than using the normal distribution 
to approximate the distribution of aggregate demands. In 
the case of three demand sources, the benefit gained from 
computing the exact distribution for the sun of demands was 
significantly lower . due to the increased computational 
effort required. However, if the numbers of repair cycles 
are small, that is not much of a problem and it is recom- 
mended that in the cases of three or four demand sources 
that the exact distribution for the sum of demands be 
computed. Otherwise for the cases of three or more demand 
sources the normal approximation can successfully be applied 
if certain parameters are held within a limited range of 
values. 

All simulation results presented in Chapter VI indicated 
that the signs of the actual deviations caused by approxi- 
mating the distribution of the aggregate demands were inde- 
pendent from the signs of the deviations resulting from 
applying the monovariable approximation formula. If, for 
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one specific parameter setting, the signs are different then 
their effects on the overall deviation get partially compen- 
sated. Otherwise, the overall deviation is about egual to 
the sum of the two single deviations. 

B. FACTORS TD BE CONSIDERED FOR PRACTICAL APPLICATIONS 

In case of binomially distributed demand sources the 
following factors should be considered. 

The NARF model incorporates mutually independent binomi- 
ally distributed demands. As was indicated in Chapter II, 
the Poisson distribution can be used to approximate the 
binomial distribution. The approximation provides an accep- 
table level of accuracy if n>20 and p<0.3 5; or if 
n>100 and np<10 [Ref. 10 p.79]. If the actual demand 
distributions at the NARF all exhibit such parameter values 
then the monovariable formula can be used along with the 
Poisson distribution for the aggregate demands. 

As a last resort, simulation of the actual demand 
behavior may be the only approach for determining the 
expected shortage costs per period. In particular, if one 
of the restrictions listed above is violated, simulation 
should be used to estimate the expected shortage costs. 

C. SUGGESTIONS FOR FURTHER ANALYSIS 

The approximation for the expected shortage costs 
appeared to be more accurate when the timing of demand moved 
to the end of the repair cycle and when the priority disci- 
pline was based on expected demand. In the case of the 
latter result, a minor change in the policy of running the 
system is suggested. 

As the stock gets more and more depleted during the 
period a rational policy would be to reserve stock for 
customers with the higher unit shortage costs and expected 
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higher demands. The higher the variability of the customer 
specific unit shortage costs and the lower the holding cost 
per stocked item, the greater is the benefit from switching 
from a "first-come, first-served" rule to a more sophisti- 
cated rule which allows for the different economic impor- 
tance of the different demand sources. Analysis has to show 
whether the costs for implementing this flexible policy are 
justified by cost reductions in the sum of economic costs 
associated with running the inventory system. 

In the case of the former result, further study of the 
influence of timing is warranted. Perhaps a shortage cost 
approximation can be developed which includes the timing as 
a parameter. 
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APPENDIX A 

COMPUTER PROGRAM FOR SIMULATING TWO BINOMIAL DEMAND SOURCES 



$J03 NAME: WILHELM F. NOEGGERATH 

DATE: JULY, 18TH 1S84 

BINOMIAL DISTRIBUTION. TWO DEMAND SOURCES 

- SIMULATION OF BERNOULLI TRIALS PER REPAIR CYCLE 

- MONOVARIABLE APPROXIMATING FORMULA 

- HYPERGEC METRIC ORIENTATED APPROXIMATION FORMULA 

- CDF ORIENTATED APPROXIMATION FORMULA 



VARIABLES: 

APPROX RESULT OF MONOVARIABLE APPROXIMATING FORMULA 
AUX. . AUXILIARY VARIABLES 
AVECOS AVERAGE COST PER PERIOD 

BIAELI ELIMINATES BIAS IF TWO DEMANDS IN ONE TIME UNIT 
AND ONE UNIT OF STOCK LEFT ONLY 
COSCDF RESULT OF CDF- ORIENTATED APPROXIMATION FORMULA 
COSCRT COST FOR CURRENT SIMULATED PERIOD 
COSCUil CUMULATED COST FOR ALL SIMULATED PERIODS 
COSHYP RESULT OF HYPERGEOM. ORIENTATED APPROXIM. FORMULA 
COSSQ CUMULATED PERIOD SQUARED COSTS 

CTR COUNTS PERIODS IN WHICH SHORTAGES OCCUR DURING 
SIMULATION 

DEVCDF DEVIATION IN PERCENT OF RESULT CDF ORIENTATED 

APPROXIMATING FORMULA FROM RESULT OF SIMULATION 
DEVKYP DEV. IN PERCENT RESULT HYPGEO. APPROXIM. FORMULA 
DEV MON DEV. IN PERCENT RESULT MONOVAR. APPROXIM. FORMULA 
DMD SCALAR. CUMULATIVE DEMAND IN ONE PERIOD 
DMD 1 ARRAY OF DEMANDS WITHIN ONE PERIOD CAUSED 
BY END ITEM #1 

DMD2 ARRAY OF DEMANDS WITHIN ONE PERIOD CAUSED 
BY END ITEM #2 

INTVL1 LENGTH OF REPAIR CYCLE FOR END ITEM #1 
INTVL2 LENGTH OF REPAIR CYCLE FOR END ITEM #2 
ISEED INITIAL SEED, GETS PRINTED 

ITMQIT INDICATES FOR WHICH END ITEM IN CURRENT SIMULATED 
PERIOD LAST DEMAND SATISFIED 
MU SUM: MU 1 + MU 2 

M U 1 MEAN BINOHIAL-DIST. DEMAND CAUSED BY END ITEM #1 

MU2 MEAN BINOMIAL-DIST. DEMAND CAUSED 3Y END ITEM #2 

N NUMBER OF DISCRETE TIME UNITS PER PERIOD 

!! **** N MUST EE A MULTIPLE OF N1 AND OF N2 **** !! 

N 1 NUMBER OF REPAIR CYCLES END ITEM #1 IN ONE PERIOD 

N2 NUMBER OF REPAIR CYCLES END ITEM #2 IN ONE PERIOD 

PI1 SHORTAGE COST ITEM #1 
PI2 SHORTAGE COST ITEM #2 

FMF PROB. MASS FCT FOE DMD PER PERIOD CAUSED 
BY E-ITEMS #1 AND #2 

PMF1 PROBABILITY MASS FCT FOR DMD PER PERIOD 
CAUSED BY E-ITEM #1 

PMF2 PROBABILITY MASS FCT FOR DMD PER PERIOD 
rAflRFR RY F-TTFM #2 

PI PROB A B. OF DEMAND CAUSED BY END ITEM #1 

IN ONE REPAIR CYCLE 

P2 PROBAB. OF DEMAND CAUSED BY END ITEM #2 

IN ONE REPAIR CYCLE 
REPLI NUMBER OF SIMULATED PERIODS 
SATDMD SATISFIED DEMAND 

SDCOS ESTIMATED STANDARD DEVIATION OF AVERAGE COST 
PER PERIOD 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



SDMEAN 

SEED 

SHEXAC 

SHSIMU 

SUM 

U 

Y 



SDCOS AS A FRACTION OF THE MEAN ESTIMATE 
SEED FOR RANDOM NUMBER GENERATION 
EXACT FRACTION OF PERIODS IN WHICH SHORTAGES 
SHOULD OCCUR 

FRACTION OF PERIODS IN WHICH SHORTAGES 
OCCURED DURING SIMULAT 

ACCUMULATES TO EXPECTED NUMBER OF STOCKOUTS 
ONE RANDOM SAMPLE VALUE FROM UNIFORM (0,1) 
INITIAL QUANTITY AVAILABLE TO SATISFY DEMAND 



DECLARATIONS * 

INTEGER REPLI, N, SEED, ISEED . INTVL 1 , INTVL2 ,S ATD MD, Y, 

* N1,N2, DMD, DMD 1(1000) ,D HD 2(1000) , ITMQIT , CTR 
REAL U (1) ,P1,P2,P±1,PI2,C0SCRT, COSCUM, COSSQ,AVECOS, 

* SDCOS, APPROX, MU, MU 1, MU 2, BIAELI, SUM, SHSIMU, 

* SDMEAN, DEVCDF,DEVMON,DEVHYP 

REAL*8 PMF (1000) , PMF1 ( 100 ),PMF2(100 ) , COS HYP , COSCDF, 

* SHEXAC 
INITIALIZATION: 



N 




20 


N 1 


= 


4 


N2 


= 


5 


PI 




0. 25 


P2 


= 


0. 75 


PI1 


= 


100. 0 


PI2 




1 000. 0 


Y 


= 


7 


REPLI 


= 


1 00 


SEED 


= 


150449 


INTVL 1 


= 


N/N 1 


INTVL2 


= 


N/N2 


ISEED 


- 


SEED 


COSCUM 


— 


0.0 


COSSQ 
CTR = 


o' 


0.0 



EACH EXECUTION OF THE FOLLOWING DO-LOOP GENERATES ONE 
VALUE FOE THE ’SHORTAGE COST PER PERIOD’ BY SIMULATING 
ONE BERNOULLI-TRIAL FOR EACH REPAIR-CYCLE 
(I.E. INTVL1, INTVL2 ) OF THE END ITEMS 

DO 300 1=1 , REPLI 
J = 1 

WHILE ( J. LE.N ) 

DMD 1 (J) = 0 

DMD2 (J) = 0 

J = J + 1 

ENDWHILE 
DMD = 0 

DO 100 J=1, N, INTVL 1 

CALL SRND (SEED , 0, 1 , 2, 0) 

IF (U (1) .LE. PI) THEN 
DMD 1 ( J) = 1 
DMD = DMD + 1 
ENDIF 

100 CONTINUE 

DO 200 J=1, N,INTVL2 

CALL SRND (SEED, U, 1 ,2, 0) 

IF (U (1) .LE. P2) THEN 
BMD2(J) = 1 
DMD = DMD + 1 
ENDIF 

200 CONTINUE 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



THE PROGRAM ASSUMES THAT THE DEMAND OCCURS RIGHT AT THE 
BEGINNING OP THE REPAIR CYCLE. THE VARIABLE ' B IA ELI ' 
COMPENSATES THE RESULTING BIAS. THE BIAS OCCURS IN 
SITUATIONS IN WHICH THERE IS ONLY ONE ITEM LEFT TO 
SATISFY DEMAND AND THERE ARE TWO DEMANDS OCCURRING 
IN THAT TIME UNIT. IN THE PROGRAM ALWAYS THE SAME 
DEMAND SOURCE GETS SERVED FIRST. 

THE ASSUMPTION THAT THE DEMAND OCCURS RIGHT AT THE BE- 
GINNING OF THE REPAIR CYCLE IS OF MORE IMPORTANCE FOR 
THE NUMERICAL RESULT OF THE SIMULATION THE SMALLER AND 
THE MORE DIFFERENT THE NUMBERS OF REPAIR CYCLES ARE. 



* 



BIAELI =0.0 
IF ((DMD.GT.Y) 
SATDMD = 0 
J = 1 
LOOP 
IF 



• AND. (Y.GT.O) ) THEN 



{ DMD1 (J) . EQ. 1 ) THEN 
SATDMD = SATDMD «• 1 
DM D1 ( J) = 0 
;IT = 1 



QUIT 

THEN 

+ 1 



ITMQ 
ENDIF 

IF ( SATDMD. EQ.Y ) 

IF j DMD2 (J) . EQ. 1 ) 
SATDMD = SATDMD 
DMD2 ( J) = 0 
ITMQIT = 2 
ENDIF 

IF { SATDMD. EQ. Y ) 

J = J + 1 
ENDLOOP 

IF ( ( DMD2 (J) . EQ. 1 ) . AND. 

BIAELI = (PI 1 + PI 2) /2 
ENDIF 
ENDIF 

COSCRT = BIAELI 



QUIT 



(ITMQIT.: 
.0 - PI 2 



!Q. 1) > 



IF 



( DM D . G T.Y 
CTR = CTE + 



:tr = 

J = 1 
WHILE fJ.LE. 
COSCRT = 

J = J 
ENDWHILE 
ENDIF 



) THEN 

1 






1 



SCRT + FLOAT (DMDIjJ) ) *PI 1 



+ FLOAT (DMD2 (J) ) *PII 



THEN 



’COSCUM’ AND ' COSSQ ' CUMULATE THE VALUES WHICH ARE LATER 
ON NEEDED FOR COMPUTATION OF THE SAMPLE MEAN AND 
ESTIMATED STANDARD DEVIATION OF THE SAMPLE MEAN. 



COSCUM = COSCUM +COSCRT 
COSSQ = COSSQ + COSCRT**2 
300 CONTINUE 
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60 1 



602 



II U 1 = P 1 AFLOAT (N1) 

MU2 = P2*FL0AT (N2 
MU = MU 1 + MU2 

COMPUTATION OF SAMPLE MEAN AS MEAN ESTIMATE AND 
COMPUTATION OF ESTIMATED STANDARD DEVIATION OF 
MEAN ESTIMATE 

AVECOS = COSCUM/FLOAT (EEPLI) 

AUX = C3 SCUM** 2/FLOAT (REPLI) 

A UX1 = FLOAT (REPLI) *FLOAT (REPLI- 1) 

SDCOS=SORT ((COSSQ-AUX) /AUX1) 

^KRITE (9,601) REPLI^N, IS EED, N 1 , P 1 , MU 1 , P1 1 , Y,N2,P2, 

^FORMATj’ ' t 7 1 ( '_*’J / t 1. REPLICA TIONS: ' £ I6 <.124, 

* • pii=* ( 

* F7.2.T67, f I=» 14,/, ' # OF REPCYCLS 2:', 13 

• P2 = ' F6.4, f MU2=',F7.2,' PI2=',F7. 



'EM AT (' ',71 '*'),/ ' REPLICATIONS: '16 ,T24, 

'TIME UNITS PER PERIOD:' ,I4,T53, ' INIT-SEED=' ,110,//, 
' # OF REPCYCLS 1:',I3,' P1=',F6.4,' MU1=',F7.2, 



WRITE (9,662)." AVECOS, SDCOs" 



2 ,/ ) 



* 

* 

* 



format}' ' ' ^feESULT"o£ BERNOULLI TRIALS FOR EACH' 



ESTIMATED 



REPAIR CYCLE: • /, 

ESTIMATED MEAN OF SHORTAGE COSTS PER', 
PERIOD' , F 8. 3, /, 

ST-DEV. FOR ESTIMATE OF THE MEAN : ' , F8. 3) 
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c 

c 

c 



604 



C 

c 



APPLYING THE MONOVARIABLE APPROXIMATING FORMULA AND 
EVALUATING THE EXACT EXPECTED FRACTION OF PERIODS 
IN WHICH SHORTAGES OCCUR 
1 = 0 

WHILE { I. LE. 



* 

* 



PHFt (1+1) 



ni ) 

=DGAMMA (DFLOAT (NI + 1) ) / 
(DGAMMA (DFLOAT (N 1-1+ 1] 



I = 1+1 
ENDWEILE 
1=0 

WHILE ( I. LE. N 2 ) 



DGAMMA (DFLOAT (1+ 1) ) ) 

PMF1 (1+1) = PMF1 (I+1)*(DBLE (Pi) ) **I - 
(DBLE (1.0 - PI) ) ** (N 1 -I) 



PMF2 (1+ 1) =DGAMM A (DF LO AT (N2 + 1 ) ) / 
'•" 2 - 1 + 1 ) ) 



) ) **I 
(N2-I) 



7 DGAMMA (DFLOAT (N2- 

* DGAMMA (DFLOAT (1+ 1) ) 

PMF2 (1+ 1) = PMF2 (1+1) *(D3LE (Pi 

* (DBLE ( 1 . 0 - P 2 ) ) **’ 

I = 1+1 

ENDWHILE 
1 = 0 

WHILE ( I. LE. N 1 +N2 ) 

PMF(I + 1) = 0.0D0 
I = I +1 
ENDWHILE 
1=0 
WHILE ( I 
J = 0 
WHILE 
PMF 
J = 

ENDWHILE 
1 = 1+1 
ENDWHILE 
SUM = 0.0 
SHEXAC = 0.0D0 
I = Y + 1 
WHILE (I.LS. NI 
SUM = SUM + 

SHEXAC = SHEXAC 
I = 1+1 
ENDWHILE 

APPROX= ( (P1 1 * MU 1 + PI 2*M U2) / ( MU 1 + MU 2) ) *SUH 
WRITE (9.604) APPROX 

FORMAT (* RESULT OF MONOVAR. 

* 'Approximating formula: • ,F8. 3,/) 



.LE.N1 ) 

[ J . LE. N2 ) 

[1+ J+ 1 ) = PMF (I+J+1) 
J + 1 



+ PM F 1 (1+1) * PMF2 (J+1) 



+ N2) 

FLOAT (I — Y) *SN3L (PMF (I 
PMF (1+1) 



♦ 1 )) 
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605 



606 



CALL HYPGEO (N 1 , N2, PMF 1 , PMF2 , P1 1 , PI2 , Y , COS H Y P) 
CALL SUBCDF (N 1 , N2, PMF 1 , PMF2 , P 1 1 , ?I2 , Y , COSCDF) 



COMPUTATION OF SUMMARY STATISTICS 



SHEXAC = SHEXAC*100.0 



SHSIMU = l 
IF (AVECOS 
DEVCDF= 
DEVKYP= 
DEVM0N= 
SDMEAN 
WRITE ( 

* 

FORMAT 



CTR/ FLOAT (RE PL I) ) *10 0.0 
.GT. 0.005) THEN 

COSCDF /AVECOS - 1.0) *100.0 
CCSHYP /AVECOS - 1.0) *100.0 
APPROX/AVECOS - 1.0) *100.0 
' (SDCOS/AVECOS) *100.0 
9,605) MU, SHEXAC, SHSIMU , SDMEAN, DEVMON, 
DEVHYP,DEVCDF 

(» SUMMARY: M0 = ’ t F5.2.» SHEXAC=' 

F6.2,'% SHSIMU=',F6. 2, ’£',//, 

• SD/MEAN=' ,F6. 1, ' % DEVMON = < , F6 . 1 , » %' 



v/ nijAW— # r q • i # / o l/£» v nun- # 

T39,*DEVHYP =, ,F6. 1,153, ' % DEVCDF=» 

EL SE 

WRITE ( 

FORMAT 



,F6. 1, " 



) 



* 

* 

* 



ENDIF 

STOP 

END 



9,606) MU, SHEXAC, SHSIMU 

(' SUMMARY: MU=',F5.2 ' SHEXAC=', 

F6.2,'5k SHSIM U=,F6.2, , 5 , / // , 

• SD/MEAN= ******* DEVMON= ********* ,T39, 
' DEVHYP= ********* ,T57, 'DEVCDF= *********) 
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c 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



601 



c 



SUBROUTINE HYPGEO (XI M AX , X2 MAX , PMF 1 , PM F2 , PI 1 , PI2 , 

* Y, COST) 

THIS SUBROUTINE HAS AS INPUT PARAMETERS: 

UMAX MAXIMAL VALUE DEMAND SOURCE#1 
X2MAX MAXIMAL VALUE DEMAND SOURCE#2 
EMF1 PR03AB. MASS ECI FOR DEMAND SOURCE# 1 , 

ALWAYS HAS TO START AT ZERO (1.3. 0,1,2, ...) 
PMF2 PROBAB. MASS FCT FOR DEMAND SOURCE#2 , 

ALWAYS HAS TO START AT ZERO (I.E. 0,1,2, ...) 

PI1 SHORTAGE COST ITEM #1 
PI2 SHORTAGE COST ITEM #2 
Y PURCHASED QUANTITY 

COST RETURNS VALUE OF HYPERGEOM. APPROXIMATION TO 
MAIN PROGRAM 

HYPGEO EVALUATES THE EXPECTED VALUE OF SHORTAGE 
COSTS ACCORDING TO THE FORMULA FOR THE EXPECTED VALUE 
OF THE HYPERGEOMETRIC DISTRIBUTION FOR EACH POSSIBLE 
COMBINATION OF XI AND X2 



THE SUBROUTINE OUTPUTS AN APPROXIMATE VALUE FOR THE 
EXPECTED SHORTAGE COST PER PERIOD 



INTEGER X1MAX, X2MAX, SHORT, Y 
E E AXi P1 1 pi2 

REAL* 8 p£lFl(100 ),PMF2(100 ),COST 
COST = O.ODO 
II = 0 



* 

* 

* 



* 



WHILE (II. LE. X1MAX) 

12 = 0 

WHILE (12. LE. X2MAX) 

IF (11+12 .GT. Y)_ THEN 
SHORT = II + 12 - Y 

COST = COST + PMF1 (11 + 1) *PMF2(I2+1) 
(DFLOAT (S HORT*1 1 ) *DBLE (PI 1) 
DFLOAT (SHORI*I 2) *DBLE (PI 2) ) / 
DFLOAT (11+12) 

ENDIF 

12 = 12 + 1 
ENDWHILE 
11 = 11+1 
ENDWHILE 

WRITE (9.601) COST 

FORMAT {’ RESULT OF HYPERGEO. APPROXIMATING 
' FORMULA: » ,F8. 3,/) 

RETURN 

END 



* 



/ 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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200 



SUBROUTINE SUBCDF ( XI H A X , X2 MA X , PHF 1 , P HF2 , P1 1 , PI2 , 

* Y, COS CDF) 

THIS SUBROUTINE HAS AS INPUT PARAMETERS: 

X 1 MAX MAXIMAL VALUE DEMAND SOURCES! 

X2MAX MAXIMAL VALUE DEMAND SOURCE#2 

PMF 1 PRCBAB.MASS FCT FOR DEMAND SOURCE# 1 , 

ALWAYS HAS TO START AT ZERO (I.E. 0,1,2, ) 

PMF2 PROB AB . MASS FCT FOR DEMAND SOURCES2 , 

ALWAYS HAS TO START AT ZERO (I.E. 0,1,2, ...) 
PI1 SHORTAGE COST ITEM #1 

PI2 SHORTAGE COST ITEM #2 

SHORT SHORTAGE #OF ITEMS 

Y PURCHASED QUANTITY 

COSCDF RETURNS VALUE OF CDF ORIENTATED APPROXIMATION 
TO MAIN PROGRAM 

THE SUBROUTINE OUTPUTS AN APPROXIMATE VALUE FOR THE 
EXPECTED SHORTAGE COST PER PERIOD 



SUBROUTINE INTERN VARIABLES: 

CDF1 CUMULATIVE PROBAB. DISTRIBUTION DEMAND SOURC 

CDF2 CUMULATIVE PROBAB. DISTRIBUTION DEMAND SCURC 

Y 1 ASSUMED SATISFIED DMD FOR END ITEM *1, 

TRUNCATED TC INTEGER 

Y 2 ASSUMED SATISFIED DMD FOR END ITEM #2, 

TRUNCATED TO INTEGER 
W1 WEIGHT FOR SHORTAGE COST PI1 

W2 WEIGHT FOR SHORTAGE COST PI2 



INTEGER XI MAX, X2MAX, SHORT, Y, Y 1 , Y2 
REAL PIl PI2 

REAL* 8 P&F 1 (100 ),PMF2(100 ),CDF1(100 ),CDF2(100 ), 
* COSCDF, W1,W2 

CDF1 (0+ 1 ) = PMF1 (0 + 1) 

DO 100 1=1 , X 1 M A X 

CDFljl+1) = CDF 1 (I) 

CONTINUE 

CDF2J0+1) = PMF2 (0+ 1) 

DO 200 1= 1 , X2M AX 

CDF2 ( 1+ 1) = CDF2 (I) 

CONTINUE 
COSCDF = 0.0D0 
II = 0 

WHILE (II. LE. X IMA X) 

12 = 0 
WHILE 



+ PMF 1 (1 + 1) 



+ PMF2 (1+1) 



IF 



(12. LE. X2MAX) 
(11+12 -GT. Y) THEN 
SHORT = II + 12 - Y 



* 

* 



ENDIF 
12 = 12 
ENDWHILE 
11 = 11 + 1 
ENDWHILE 



Y 1 =1 NT ((FLOAT (II) /FLOAT (11+12) ) * FLO A T ( Y) ) 
W1=( (CDF 1 (I1 + T)-CDF1 { Y 1 + 1 ) ) /CDF 1 (11+1) ) 

* DFLOAT (I 1 ) 

Y2 = I NT (?FLOAT (12) /FLOAT ( I1 +12) ) * FLOAT (Y) ) 
W 2= ( (CDF2 (I2 + T) -CDF2 (Y2+ 1) ) /CDF2 (12+ 1) ) 

*DFLOAT(I2) 

COSCDF = COSCDF + PMF 1 (1 1 + 1) *PMF2 (1 2+ 1) * 

DFLOAT (SHORT) * 

( (W1*PI1 + W2*PI2) / (W1 + K2) ) 



+ 1 



601 

* 



SENTRY 

SEND 



WRITE 

FORMA' 

RETURN 

END 



(9 <‘ 



6 01) COSCDF 

RESULT OF 
FORMULA: • ,F8. 3,/) 



CDF APPROXIMATING 1 , 
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appendix B 

COMPOTES PROGSAM FOR SIMULATING N BINOMIAL DEMAND SOURCES 



$JOB 



NAME: WILHELM F. NOEGGERATH 
DATE: JULY, 21ST 1984 



BINOMIAL DISTRIBUTION, N DEMAND SOURCES 

- SIMULATION OF BERNOULLI TRIALS PER REPAIR CYCLE 

- APPLYING MONOVARIABLE APPROXIMATION FORMULA 

BY USING NORMAL APPROXIMATION TO 3IN0MIAL DIST* S 



V ART ART FR • 

. AUX. AUXILIARY VARIABLES 

AVECOS AVERAGE COST PER PERIOD 

COSCRT COST FOR CURRENT SIMULATED PERIOD 

COSCUM CUMULATED COST FOR ALL SIMULATED PERIODS 

COSSQ CUMULATED PERIOD SQUARED COSTS 

CRTDMD VECTOR WHICH STORES ORIGINAL DEMAND OF CURRENT 
TIME UNIT 

CUMDMD CUMULATES DEMAND DURING PERIOD 

CTR COUNTS PERIODS IN WHICH SHORTAGES OCCUR DURING 
SIMULATION 

DEVMON DEVIATION OF RESULT OF MONOVARIABLE APPROXIMATION 
FOR EXPECTED SHORTAGE COST PER PERIOD FROM 
APPROXIMATION RESULT IN PERCENT 
DEVNSH DEVIATION OF RESULT OF MO NOV ARI ABLE APPROXIMATION 
FOR EXPECTED NUMBER OF SHORTAGES PER PERIOD FROM 
SIMULATION RESULT IN PERCENT 
DEVPIA DEVIATION OF THE AVERAGE OF THE UNIT SHORTAGE 
COST AS USED IN THE NORMAL APPROXIMATION FROM 
THE AVERAGE UNIT SHORTAGE COST AS OBSERVED 
DURING SIMULATION 

DMD DEMAND MATRIX FOR ONE PERIOD 

(ROW=DMD-SOURCE, COL=T IHE-U NI T) 

INTVL VECTOR WITH LENGTH OF REPAIR CYCLES FOR 
DIFFERENT DMD SOURCES 

N NUMBER OF DISCRETE TIME UNITS PER PERIOD 

*** N MUST BE A MULTIPLE OF EACH ELEMENT OF REPCYC *** 
NOFDI1D NUMBER OF DIFFERENT DEMAND SOURCES 
NOFSHN AVERAGE NUMBER OF SHORTAGES PER PERIOD 
DERIVED FROM NORMAL APPROXIMATION 
NOFSHS AVERAGE NUMBER OF SHORTAGES PER PERIOD 
OBSERVED DURING SIMULATION 
P VECTOR , PROBABILITIES OF DEMAND 

IN ONE REPAIR CYCLE 

PI VECTOR , UNIT SHORTAGE COSTS FOR DIFFERENT 

DEMAND SOURCES 

PINORA MEAN WEIGHTED AVERAGE UNIT SHORTAGE COST 
PISIMA AVERAGE UNIT SHORTAGE COST AS OBSERVED DURING 
SIMULATION 

REPCYC VECTOR WITH #0F REPAIR CYCLES PER PERIOD 
REPLI NUMBER OF SIMULATED PERIODS 

SATDMD CUMULATES SATISFIED DEMAND FOR CURRENT PERIOD 
SDAVE ESTIMATED STANDARD DEVIATION FOR THE MEAN 

ESTIMATE FOR THE EXPECTED SHORTAGE COST PER 
PERIOD DERIVED FROM SIMULATION DIVIDED BY THE 
ESTIMATE I.E. DIVIDED BY THE SAMPLE MEAN 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



SDCOS ESTIMATED STANDARD DEVIATION OF SIMULATED 
AVERAGE SHORTAGE COST PER PERIOD 
SDNOSH ESTIMATED STAND. - DE VIATI ON OF 
NUMBER SHORTAGES PER PERIOD 
SEED SEED FOR RANDOM NUMBER GENERATION 
SHORT #OF SHORTAGES FOR CURRENT SIMULATED PERIOD 
SHCUM CUMULATED #OF SHORTAGES FOR ALL SIMULATED PERIODS 
SHSIMU FRACTION OF PERIODS IN WHICH SHORTAGES 
OCCURED DURING SIMULATION 
SHSQ CUMULATED PERIOD SQUARED #OF SHORTAGES 
U ONE RANDOM SAMPLE VALUE FROM UNIFORM (0,1) 

UNTDMD CUMULATES DEMAND FOR ONE TIME UNIT 
Y INITIAL QUANTITY AVAILABLE TO SATISFY DEMAND 



DECLARATIONS 
INTEGER CTR 

* 

REAL 

* 

* 

* 



C 

C 

c 

c 

c 

c 

c 



TR.CRTDMD (50) ,INTVL(50) , N , NOFDM D, R EPC YC (50) , 
REPLI.SATDMD, SEED , UNTDMD, UNTSAT, XM AX , Y , CUi-lD MD 
AVECOS , COS CRT, COSCUM, COSSQ.DMD (5 0, 100) ,P (50) , 
PI (50) , SDCOS, SDNOSH, 5HCRT, SHCUM, SHSIMU, SHSQ, 
MU, U ( 1) ,NOFSHN,NOFSHS, APPROX, PI NOR A, S DAVE, 
PISIMA,CEVMON, DEVNSH, DEVPIA 



WITH THE CURRENT DIMENSION THE PROGRAM IS ABLE TO HANDLE 
50 DEMAND DISTRIBUTIONS FOR A PERIOD OF 100 TIME UNITS. 
FOR SETTINGS OF MORE THAN 50 DISTRIBUTIONS DR MORE THAN 
100 TIME UNITS THE DIMENSION STATEMENTS HAVE TO BE 
CHANGED 



SEED 





READ (5, 


10 


FORMAT ( 




1=1 




WHILE (I 




READ 


20 


FORMA 




1 = 1 + 1 




ENDWHILE 



255 

NOFDMD, 



N, Y , R EPLI 



E. NOFDMD1 
,20) REPCYCfl) 
(2 1X,I3,9X, F6. 






UNITS= ' , 13, 
! // 1 



:) \ *pm 

, :)) * (l.o-P(i) ) 

I,REPC^C (I) , I, P (I) , I , AUX, I , PI (I) , 



WRITE (9.30) NOFDMD, N,F.EPLI, SEED 
30 FORMAT (* • , • N OFDMD= • , 12 , 1 #OF TIME 

* • R EP LI= ’ , 15 , • I-SEED = ' ,18.,/ , 

* T2,'#OF REP-' ,T1 3, ' PROBA-* ,T38.‘ UNIT-SHORT. ' , 

* T55, 'FIT' , /. 12 , 'CYCLES' ,T13, ' BxLIT Y ' ,T25, 

* 'MEAN' , T38 , 'COST' , T 53 , ' N*P' , T5 8, • N*Q ' ,/) 

WHILE (I . LE. NOFDMD) 

AUX = FLOAT (REPCYC (I) 

AUX1 = FLOAT (REPCYC (I 
WRITE (9,40) I, REPCYC 

* AU X A U X 1 

40 FORMAT (• ' , 12, ' : • , 13 ,16 , • : ',F4.2,I6 ': • 

* F6.2.I5,': ' , F7 . 2, 3X , f 1 5. *1 , F5 . S ) 
INTVL(I) = N/REPCYC (I) 

1=1+1 

ENDWHILE 

SUMMU = 0.0 

SUMVAR = 0.0 

VM JY = n 

DO 50 1=1, NOFDMD 

MU = FLOAT (REPCYC (I) ) *P (I) 

SUMMU = SUMMU + MU 
SUMVAR = SUMVAR + MU* ( 1 . 0-P (I) ) 

XMAX = XMAX + REPCYC (I) 

50 CONTINUE 

STDDEV = SORT (SUMVAR1 

WRITE (9,60) SUMMU, STDDEV, Y, XMAX 
60 FORMAT (' ', /.'TOTAL MEA N= ' , F6 . 2, ' 



F5-2, 



Y= ' ,13, 



TOTAL S-DEV= ' , 



MAX- DM D=' , 14) 
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COSCUM = 0.0 
COSSQ = 0.0 
Sh’CUM = 0.0 
SHSQ = 0.0 
CTR = 0 
c 

C EACH EXECUTION OF THE FOLLOWING DO-LOOP GENERATES ONE 
C VALUE FOR THE 'SHORTAGE COST PEE PERIOD' 31 SIMULATING 

C ONE BERNOULLI-TRIAL FOR EACH REPAIR-CYCLE (I. E. INT VL (I) ) 

C OT THE END-ITEMS 
C 

DO 900 1 = 1 , REP LI 

DO 200 J= 1 . NOFDMD 
DO 100 K = 1,N 

DMD (J ,K) = 0.0 
100 CONTINUE 

200 CONTINUE 

CUM DMD = 0 
DO 400 J=1, NOFDMD 
IAUX = INTVL(J) 

DO 300 K= 1 , N, IAUX 

CALL SEND (SEED , U , 1 , 2 , 0) 

IF ( U (1) . LE. P(J) ) THEN 
DMDjJ.K) = T.O 
CUMDMD = CUMDMD + 1 
ENDIF 

300 CONTINUE 

400 CONTINUE 

IF ((CUMDMD. GT. Y) . AND. (Y.GT.0) ) THEN 
SATDMD = 0 
K= 1 
LOOP 

UNTDMD = 0 
DO 500 L=1, NOFDMD 

IF ( DMD(L,K) .EQ. 1.0 ) THEN 
UNTDMD = UNTDMD + 1 
CRTDMD(L) = 1 

FT 

CRTDMD(L) = 0 
ENDIF 

500 CONTINUE 

DO 600 J= 1 , NOFD MD 

IF (DMD (J,K) .EQ. 1.0) THEN 
SATDMD = SATDMD + 1 
DMD ( J, K) = 0.0 
ENDIF 

600 CONTINUE 

IF (SATDMD. GE. Y) QUIT 
K = K + 1 
ENDLOOP 



IF MORE THAN THE AVAILABLE QUANTITY GOT ASSIGNED 
TO THE CUSTOMERS A CORRECTION HAS TO BE PERFORMED. 



ADDITIONALLY THE EIAS CAUSED BY CHECKING THE DEMAND 
MATRIX ALWAYS IN THE SAME SEQUENCE HAS TO BE ELIMINATED. 

IF fSATDMD. GT. Y) THEN 
DO 750 J=1. NOFDMD 

IF (CRTDMD (J) . EQ. 1) THEN 

DMD ( J, K) = FLOAT (SATDMD - Y) / 

* FLOAT (UNTDMD) 

ENDIF 

750 CONTINUE 

ENDIF 
ENDIF 



73 



800 

850 



900 

C 

C 



950 



C 



C 

c c 
c 



960 



970 



C 



COSCRT = 0.0 
SKCRT =0.0 

IF ( CUMDMD .GT. Y ) THEN 
CIS = CTR + 1 
DO 850 J= 1, NOFDMD 
DO 800 K=1,N 

COSCRT = COSCRT + DMD ( J, K) *PI ( J) 
SHCRT = SHCRT + DMD(J,K) 

CONTINUE 

CONTINUE 

ENDIF 

COSCUM = COSCUM +COSCRT 
COSS2 = COSSQ + COSCRT**2 
SHCUM = SHCUM + SHCRT 
SHSQ = SHSQ + SHCRT** 2 
CONTINUE 



AVECOS = COSCUM/FLOAT (REPLI) 

AUX2 = COSCUM**2/FLOAT (REPLi) 

AUX3 = FLOAT (REPLI) *FLOAT (REPLI-1) 

SDCOS=SQRT ( (COSSQ-AUX2) /AUX3) 

NOFSHS = SHCUM/FLOAT (REPLI) 

A UX4 = SHCUM**2/FLOAT (REPLI) 

SDNOSH=SQRT ( (S HSQ- AUX4 )/AUX3) 

SHSIMU= (FLOAT (CTR) /FLOAT(REP) 

WRITE (9,950) AVECOS . NOFS HS , SL. 

FORMAT{’ SIMULATION: AVE-COST= ’ , F8. 2, 

* ’ #OF SHOR”S=’,F9.4,/, 'SHSIMU = ' ,F6.2, 

* '* SD-COST= ' , F8 . 2 , ' SD-#OF SH ’ ’ S=' , F8 . 4) 



‘LI) ) *100. 0 
NOFS fiS,SHSIMU,SDCOS, SDN OSH 



CALL NORMAL (Y, NOFDMD, R EPC YC , P , PI , NCF SH N, ST DD EV , SU MM U, 
* XMAX,PINORA, APPROX) 

OMPUTATION OF SUMMARY STATISTICS 



IF (NOFSHS. GT. 
PISIMA = AV 
IF ((AVECOS 
SDAVE = 
DEV MON = 
DEVNSH = 
DEVPIA = 
WRITE (9 
s= 

FORMAT (’ 



* 

* 

* 



ELSE 

WRITE 
ENDIF 
ELSE 

WRITE 
FORMAT ' ( 

* 

ENDIF 

STOP 

END 



0.00005) THEN 

ECOS / NOFSHS 

.GT.0.0 0 5) .AND. (PISIMA. GT. 0 . 005) ) THEN 
SDCOS / AVECOS) *100.0 
APPROX/A VECCS - 1.0) *130.0 
NOFSHN/NOFSHS - 1.0) *100.0 
PINORA/PISIMA - 1.0) *100.0 
60) SDAVE, PISIMA, PINORA, DEVMON ,DEVN SH , 
DEVPIA 

’/, ’SUMMARY STATISTICS:’,/, 

T1 , ' SD AVE= ’ ,F9 . 2 , ’ , T20 , ’ PISIMA= ’, 

F8. 2, ’ $’ ,T40, ’ PINORA =' ,F8. 2,' $ ’ ,/,T1, 
’DEVMON=',F8.2, , %' ,T20, * DEVNSH = ’ , F3 . 2 , 
’S»,T40,»DEVPIA=*,F8.2,’S* ) 

(9,970) 



( 9 , 970 ) 



,/,’NO SURVEY STATISTICS 
/, ’AT LEAST ONE DIVISOR 



COMPUTED, ' , 
TOO SMALL’) 
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c 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



c 



SUBROUTINE NORMAL (Y, NOFDMD, RZPCYC,? , PI, ND FSHN , SIDDEV , 
* 5 U K M U , X M A X , ? I N 0 R A , A 2 ? R 0 X ) 

THE SUBROUTINE 'NORMAL* COMPUTES AN EXPECTED VALUE FOR 
THE SHORTAGE COST PER PERIOD ACCORDING TO THE MGNO- 
V ARIABLE APPROXIMATION FORMULA 31 USING A NORMAL 
DISTRIBUTION AS AN APPROXIMATION FOR THE DISTRIBUTION 
OF THE SUM OF DEMANDS PER PERIOD. 



SUBROUTINE INTERN VARIABLES: 

CDFL STANDARD NORMAL DISTRIBUTION EVALUATED AT LL 

CDF F. STANDARD NORMAL DISTRIBUTION EVALUATED AT RL 

CRT CURRENT SUM OF DEMANDS ( GREATER THAN Y ) 

LL (CRT - 0.5) TRANSFORMED TO STANDARD NORMAL 

RL (CRT + 0.5) TRANSFORMED TO STANDARD NORMAL 

SHORT CURRENT NUMBER OF SHORTAGES 
SUMPI SCALAR, CUMULATED ?(I)*?I(I); 1 = 1, ...,N 

riFfT jsi^TnMq 

INTEGER NOFDMD, Y,REPCYC (9) , XMA X 

REAL APPROX, SUMMU, MU, SU MV Art, SUM PI, PI (50) ,? (50) , 

* SHORT, SIDDEV, NOFSHN, Z, PINORA, CRT, 

* CDFL,CDFR,LL, RL 

SUMPI = 0.0 
DO 1 1=1, NOFDMD 

SUMPI = SUMPI + FLOAT (REPCYC (I) ) *P (I) *PI (I) 

1 CONTINUE 

SHORT = 1.0 
CRT = FLOAT (Y+ 1 ) 

NOFSHN =0.0 

LL = (CRT -0.5 -SUMMU) /SIDDEV 

THE ISML LIBRARY SUBROUTINE ' MDNOR ( A, 3 )* EVALUATES 
THE STANDARD NORMAL DISTRIBUTION AT A AND passes IIS 
VALUE VIA THE PARAMETER 3 TO THE CALLING PROGRAM 



CALL MDNOR ( LL, CDFL) 

RHILS ( CRT. LE. FLOAT (XMAX) ) OO 

RL = (CRT + 0.5 - SUMMU) /SIDDEV 
CALL MDNOR ( RL, CDF R) 

NOFSHN = NOFSHN + SHGRI * ( CDFR - CDFL) 
CDFL = CDFF. 

SHORT = SHORT + 1.0 
CRT = CRT + 1.0 
ENDRHILE 

PINORA = SUMPI/SUMM U 
APFF.OX = PINCRA * NOFSHN 
WRITE (9.2) APPROX, NOFSHN 
2 FORMAT (» * ,/, * NOEHAL-APPROXIM. : ' , 

* ’ AVE— COST= ', F8. 2, * SOF SHOR"S=* 

RETURN 
END 

SENTRY 



F 9. 4) 
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APPENDIX C 

SIMULATION RESULTS FOB THREE DEMAND SOURCES, Y INCREASING 



NOFDM D= 3 
#OF REP- 



#0F TIME UNI TS= 72 
PROBA- 



REPLI =20000 I-SEED= 901941 
UNIT-SHORT. FIT 



CYCLES 


BILITY 


MEAN 


COST 


N*P 


N*Q 


1 : 18 


1: 0.20 


1 : 


3.60 1: 


300.00 


3. 6 


1 4.4 


2: 36 


2: 0.40 


2: 


14.40 2: 


300.00 


14.4 


2 1.6 


3: 12 


3: 0.30 


3: 


3.60 3: 


100.00 


3.5 


8.4 


TOTAL MEAN = 


21 . 60 


TOTAL S- 


DE V= 3.75 


Y= 20 


M AX-DMD= 


6 6 



SIMULATION : 
SHSIM U= 61 . 06 % 



AV E-COST= 655. 11 
SD-COST= 5.27 



#0 F SHOR* S = 
SD- #OF SH ’ 3 = 



2.4178 

0.0195 



NORMAL-APPROXIM. : AVE-COST= 646.67 #OF SHOR'S= 2.4250 



SUMMARY STATISTICS: 

SD/A VE= 0.80% PISIMA= 270. 95$ 

DE VMON= -1.29% D EVNS H= 0.30% 



PI NORA= 266.67$ 
DEVPIA= -1.58% 



NOFDMD= 3 #OF TIME UNITS= 72 REPLI=20000 I-SEED= 901941 



#OF REP- 
CYCLES 

1 : 18 
2: 36 
3: 12 



PROBA- 
BILITY 
1 : 0.20 
2: 0.40 
3: 0.30 



TOTAL MEAN= 21.60 



MEAN 
1 
2 
3 



3.60 

14.40 

3.60 



UNIT-SHORT. 
COST 



TOTAL S- DE 7= 3.75 



300.00 
300.00 

100.00 

Y= 23 



FIT 

N*P H*Q 
3.6 14.4 
14.4 21.6 
3.6 8.4 

MAX-DMD= 66 



SIMULATION: AVE-COST = 246.61 

SHSIM U= 30.36% SD-COST= 3.37 



#OF SHOR' S= 0.9046 
SD- #0F SH’S= 0.0124 



NORMAL-APPROXIM. 



A V E— COST= 238.36 #CF SHOE*S = 



0.8938 



SUMMARY STATISTICS: 

SD/AVE= 1.36% P ISIMA= 272.63$ 

DE VMO N= -3.35% D EVNSH= -1.19% 



PI NORA= 266.67$ 
DEVPIA= -2.19% 



NOFDMD= 3 #0F TIME UNITS= 72 REPLI=20000 I-SEED= 901941 

FIT 

N*P N* 
3.6 14. 
14.4 21. 

3.6 8, 

MAX-DMD= 66 



#OF 


REP- 


PROBA- 






UNIT-SHORT 


CYCLES 


BILITY 


MEAN 




COST 


1 : 


18 


1: 0.20 


1: 


3.60 


1: 300.00 


2: 


36 


2: 0.40 


2: 


14.40 


2: 300.00 


3: 


12 


3: 0.30 


3: 


3.60 


3: 100.00 


TOTAL 


MEAN = 


21.60 


TOTAL S- 


DE V= 3. 


75 Y= 26 



SIMULATION: AV E-COST= 63.14 

SHS IM U= 9.80% SD-COST= 1.66 

NORMAL-APPROXIM.: AVE-COST= 

SUMMARY STATISTICS: 

SD/A VE= 2.64% P I SI M A = 275.29$ 

DE VMO N= -7.36% D EVNSH= -4.37% 



#0 F SHOR’S= 0.2293 
SD- #OF SK'S= 0.0061 



58.49 #OF S HOR ' S= 0.2193 



PINORA= 266.67$ 
DEVPIA= -3.13% 
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APPENDIX d 

SIBULATION IESOLTS FOR THREE DEMAND SOORCES 



NOFDMD= 3 # OF TIME DNITS= 60 REPLI=20000 I-SEED= 



1809 



#OF REP- 
CYCLES 

1: 10 
2: 12 
3: 15 

TOTAL ME AN= 



PROBA- UNIT- SHORT. 

BILITY MEAN COST 

1 : 0.20 1 : 2.00 1 : 100.00 

2; 0.30 2: 3.60 2: 300.00 

3: 0.80 3: 12.00 3: 200.00 

17.60 TOTAL S-DEV= 2.55 Y= 20 



FIT 

N*P N*Q 
2.0 8.0 
3.6 3.4 

12.0 3.0 

MAX-DMD= 37 



SIMULATION. 
SHS IM U= 12.6135 



AV r-COST= 
SD-COST= 



49.11 #OF SHOE' S= 0.2333 
1.09 SD— #OF SH * S = 0.0051 



NCRMAL-APPROXIM. : AVE-COST= 



48.93 #0 F S HOR ' S = 



0.2340 



SUMMARY 
SD/A V E= 
DE VMON= 



STATISTICS, 

2 . 22 % 



P I SI M A = 
D EVNS H= 



210.47$ PINORA= 209.09$ 

0.29% DEVPIA= -0.65% 



NOFDMD= 3 #OF TIME UNITS= 72 REPLI=2Q000 I-SEED= 94009 



#OF 


EEP- 


PROBA- 




UNIT-SHORT. 


FIT 


CYCLES 


BILITY MEAN 




COST 


N*P N*Q 


1: 


12 


1: 0.40 1: 


4.80 


1: 


300.00 


4.8 7.2 


2: 


24 


2: 0.50 2: 


12.00 


2: 


300.00 


12.0 12.0 


3: 


18 


3: 0.60 3: 


10.80 


3: 


100.00 


10.8 7.2 


TOTAL 


ME AN = 


27.6 0 TOTAL S- 


DEV= 3.63 




Y= 32 


M AX-DM D= 54 




SIMULATION: AV E-COST= 


45.55 


#OF 


SHOR' 5= 


0.2010 


SHS IM U= 9. 


05% SD— COST= 


1.24 


SD- 


#0 F SH’S 


= 0.0054 


NORMAL-APPROXIM.: AVE-COST= 


43.69 


#OF 


SHOR' S= 


0. 1970 



SUMMARY STATISTICS: 

SD/A V E= 2.72% PISIMA= 226. 57$ PINORA= 221.74$ 

DEVMON= -4.08% D EVNSH= -1.99% DEVPIA= -2.13% 



NOFDM D= 3 #OF TIME UNITS= 60 REPLI=20000 I-SSED= 



891 



#OF REP- 
CYCLES 
1: 15 

2 : 20 
3: 12 

TOTAL M E AN : 



PROBA- 
BILITY 
1: 0.70 
2; 0.30 
3: 0.60 
23.70 



MEAN 
1 
2 
3 



10. 50 
6.00 
7.20 



TOTAL S-DE V= 3.20 



UNIT-SHORT. 

COST 

1 ; 100.00 

2: 300.00 

3: 300.00 

Y= 27 



FIT 

N*P N*Q 
10.5 4.5 

6.0 14.0 
7.2 4.8 

MAX-DMD= 47 



SIMULATION 
SHS IM U= 11.55% 



AV I-COST= 
SD-COST= 



NORMAL-APPROXIM. : AVE-COST= 



49.95 #OF S HOR ' S = 0.2356 

1.22 SD-IOF SH'S= 0.0055 



52.36 #OF S HOR ' S = 



0. 2477 



SUMMARY 
SD/A VE= 
DEVMCN= 



STATISTICS; 
2. 44% 
4.83% 



P ISIMA = 
DEVNSH= 



211 . 

5. 



97$ 

1 1 % 



PINORA= 

DEVPIA= 



211. 39$ 
-0.27% 
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NOFDMD= 3 


#OF TIME 


UNITS= 60 REPLI = 


; 20000 I-S E ED = 


312774 


#OF 


REP- 


PROBA- 




UNIT-SHORT. 


FIT 


CYCLES 


BILITY 


MEAN 


COST 


N*P N*Q 


1 : 


10 


1: 0.30 


1: 3.00 


1: 300.00 


3.0 7.0 


2: 


15 


2: 0.60 


2: 9.00 


2: 300.00 


9.0 6.0 


3: 


20 


3: 0.50 


3: 10.00 


3: 100.00 


10.0 10.0 


TOTAL 


MEAN 


= 22.00 


TOTAL S-DE V= 3.27 


Y= 27 MAX 


- D M D = 45 



SIMULATION: AVE-COST= 17.29 

SHS IM 0= 4.88% SD-COST= 0.68 

NORMAL-APPROXIM.: AVE-COST= 

SUMMARY STATISTICS: 

SD/A VE= 3.93% PISIMA= 191.12$ 

DEVMON= 6.71% D EVNS H= -2.46% 



#0 F SHOE 1 S = 
SD-#OF SH'S= 



18.45 #OF SHOE ’ S= 



0.0904 

0.0033 

0.0382 



PINORA= 

DEVPIA= 



209.09$ 

9.40% 



NOFDM D= 3 


#OF TIME 


UNI TS= 88 REPLI = 


■20000 I-S E 


1 II 
1 Q 
i to 


932675 


#OF REP- 


PROBA- 






UNIT-SHORT. 




FIT 


CYCLES 


BILITY 


MEAN 




COST 




N*P N*Q 


1 : 11 


1 : 0.10 


1 : 


1.10 


1: 100.00 




1.1 9.9 


2: 44 


2: 0.30 


2: 


13.20 


2: 100.00 




13. 2 30. 8 


3: 22 


3: 0.90 


3: 


19.80 


3: 300.00 




19.8 2.2 


TOTAL MEAN 


= 34. 1 0 


TOTAL S- 


DE V= 3.49 


Y= 39 


MAX 


-DMD= 77 


SIMULATION: AVI-COST= 


24.60 


If OF S HO R ’ 5 = 




0.1268 


SHSIM U= 5 


.86% 


SD-COST= 


0.85 


SD- #OF SH'S 


= 


0. 0043 


NORMAL-APPROXIM.: AVE-COST= 


27. 17 


#OF SHOR’S= 




0. 12 57 


SUMMARY STATISTICS: 












SD/A VE= 


3.47% 


PISIMA= 


193. 93$ 


PI NOR A= 


216 


. 13$ 


DEV MO N= 


10.45% 


D EVNS H= 


-0. 89% 


DEVPIA= 


1 1 


. 45% 



NOFDMD= 3 #OF TIME UNITS= 84 REPLI=20000 I-SEED= 928698 



#OF REP- 
CYCLES 

1 : 21 
2: 42 
3: 14 



PROBA- 

BILITY 

1: 0.30 
2: 0.60 
3: 0.20 



MEAN 



UNIT-SHORT. FIT 

COST N*P N*Q 



6.30 

25.20 

2.80 



TOTAL MEAN = 34.30 TOTAL S-DEV= 4.09 



300.00 

300.00 

100.00 

Y= 4 0 



6. 3 14. 7 
25.2 16.8 
2.3 1 1.2 

MAX-DMD= 77 



SIMULATION: AV E-COST= 39.55 

SHSIM U= 6.12% SD-COST= 1.31 



If OF SHOE ' S = 
SD- #OF SH'S= 



NORMAL-APPROXIM.: AVE-COST= 



42.72 JfOF SHOR , S = 



0.1366 
0.00 46 

0. 1506 



SUMMARY STATISTICS: 

SD/AVE= 3.32% PISIMA= 289.46$ PINORA= 283.67$ 

DE V MO N= 8.00% D EVNSH= 10.20% DEVPIA= -2.00% 
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NOFDMD= 


3 #OF TIME UNITS= 


48 REPLI= 


: 20000 I— S E ED = 


33 11 


# 0 F REP 


PROEA- 




UNIT-SHORT. 


FIT 


CYCLES 


BILITY MEAN 


COST 


N*P M*Q 


1: 12 


1: 0.80 1: 


9.60 


1 


100.00 


9.6 2.4 


2: 24 


2: 0.40 2: 


9.60 


2 


250.00 


9.6 14.4 


3: 16 


3: 0.10 3: 


1.60 


3 


250.00 


1.6 14.4 


TOTAL ME 


AN= 20.80 TOTAL S 


- DE V= 3.02 


X 

<13 

CM 

II 

r>4 


- DM D= 5 2 


SIMULATION: AVE-COST= 


43. 50 


#0F SHOR 'S= 


0.2196 


SHSI MU= 


10.857* S B-COS T= 


1.05 


3D- #0 F SH'S = 


0.0053 


NORMAL-APPROXIM.: AVE-COST= 


39.98 


#OF S HOR ' S = 


0.2212 


SUMMAEY 


STATISTICS: 










SD/A VE= 


2.42% PISI MA= 


198. 10$ 




P I N 0 R A = 180 


.77$ 


DEVMON= 


-8.09% DEVNSH= 


0. 72% 




DEVPI A= -8 


. 75% 




NO FD MD= 


3 #OF TIME UNITS= 


80 REPLI= 


•20000 I— SEED = 


3951 13 


# OF REP 


PROBA- 




UNIT-SHORT. 


FIT 


CYCLES 


BILITY MEAN 


COST 


N*P N*Q 


1: 20 


1: 0.90 1: 


18.00 


1 


100.00 


18.0 2.0 


2: 80 


2: 0.30 2: 


24.00 


2 


100.00 


24.0 56. 0 


3: 40 


3: 0.90 3: 


36.00 


3 


300.00 


36. 0 4.0 


TOTAL MEAN= 73.00 TOTAL S 


-DEV= 4.71 




Y = 82 MAX 


- DM D= 140 


SIMULATION: AVE-COST= 


102. 12 


#OF SHOR' S= 


0. 5370 


SHSI MU= 


17. 13% SD-COST= 


1.98 


SD-#0 F SH'S = 


0. 01 05 


NORMAL-APPROXIM. : AVE-COST= 


99.36 


#OF SHO R ’ S = 


0.5167 


SUMMARY 


STATISTICS: 










SD/A VE= 


1.94% PISIMA= 


190. 17$ 




PINOEA= 192 


.31$ 


DEVMON= 


-2.70% D EVNSH= 


-3. 78% 




DEVPI A= 1 


. 12% 



NOFD MD= 3 #OF TIME UNIIS= 60 REPLI=20000 I-SEED= 7655599 



# OF REP- 
CYCLES 

1 : 10 

2: 30 
3: 15 



PROBA- 

BILITY 



1 : 

2: 

3: 



0.20 

0.40 

0.20 



MEAN 

1 : 

-) • 
a- • 

3: 



UNIT-SHORT. 

COST 



2.00 

12.00 

3.00 



100.00 
100.00 
300. 00 



TOTAL ME AN= 17.00 



TOTAL S— DE V= 3.35 Y= 22 



FIT 

N*P N*Q 

2.0 8.0 

12.0 18.0 

3.0 12.0 

MAX-DMD= 55 



SIMULATION: AVE-COST= 
SHSI MU= 5.48% SD-COST= 

NORMAL-APPROXIM. : AVE-COST= 



14.71 #OF SHOR ' S= 0.1114 
0.55 S D- #0 F SH’S= 0.0039 



13. 23 #OF SHOR' 5= 



0.0973 



SUMMARY STATISTICS: 

SD/A VE= 3.72% PISIMA= 131.97$ 

DEVMON= -10.03% D EVN SH= -12.24% 



PINOR A= 135.29$ 
DEVPI A= 2. 52% 
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NOFDMD= 3 #0F TIME UNITS= 60 

PROBA- 

MEAN 



REPL 1=20000 I-SE2D= 



#DF REP- 
CYCLES 



PROBA- 

BILITY 



UNIT-SHORT. 

COST 



1 : 20 
2: 30 
3: 15 



1 : 0.20 
2: 0.30 
3: 0.60 



4.00 

9.00 
9.00 



300.00 

300.00 

100.00 



TOTAL ME AN= 22.00 

SIMULATION : 
SHSIMU= 10.92% 

NORMAL- APPROXIM. : 



TOTAL S-DEV= 3.62 Y= 26 



3 48 76 
FIT 

N*P N*Q 

4.0 16.0 

9.0 21.0 

9.0 6.0 



MAX-DMD= 65 



SD/AVE= 

DE7MON= 



2.49% 

- 8 . 20 % 



A V E-CO 5T= 
S D-COST= 


57. 84 

1.44 


#OF SHOR ’ S= 0.2477 

SD- #0 F SH’S= 0.0061 


A VE-COST= 


53. 10 


#OF SHOR* S= 0.2434 


’PISIMA* 
D EVNSH = 


233. 45$ 

-1. 77% 


PINORA= 21 8. 183 
DEVPIA= -6.54% 



NOFDMD= 3 #OF TIME UNITS = 72 REPLI=20000 I-SEED= 176532 



# OF REP- 
CYCLES 



PROBA- 

BILITY 



MEAN 



UNIT-SHORT. 

COST 



1: 


12 


Is 


0.90 


1: 10.80 1: 


2: 


36 


2: 


0.40 


2: 14.40 2: 


3: 


18 


3: 


0.50 


3: 9.00 3: 


)TAL 


MEAN = 


34. 


20 


TOTAL S-DEV= 3.77 



100.00 

100.00 

300.00 

Y = 39 



F IT 

N*P A N*Q 

10.8 1.2 

14.4 21.6 

9.0 9.0 

MAX-DMD= 66 



SIMULATION: AVE-COST= 26.99 
SHSIMU= 8.02% S D-COST= 0.80 



#OF SHOR 1 S= 0. 1781 
SD- #0 F SH'S= 0.0051 



NORMAL-APPROXIM. : AVE-COST= 



27.45 #OF SHOR’S= 



0. 1799 



SUMMARY STATISTICS: 

SD/A VE= 2.96% PISIMA= 151.55$ ?INORA= 152. 63$ 

DEVMON= 1.71% DEVNSH= 0.99% DEVPIA= 0.71% 



NOFDMD= 3 


#OF TIME UNITS= 78 REPLI= 


■20000 I-SEED= 3480610 


#DF REP- 


PROBA 


— 


UNI 


T-SHORT. 


FIT 


CYCLES 


BILITY MEAN 


COST 


N*P N*Q 


1: 13 


1: 0. 


40 1: 5.20 


1: 


100.00 


5.2 7.8 


2: 78 


2: 0. 


20 2: 15.60 


2: 


100.00 


15.6 62.4 


3: 39 


3; 0. 


70 3: 27.30 


3: 


300.00 


27. 3 1 1. 7 


TOTAL MEAN= 


4^ 

CD 

• 

O 


TOTAL S-DEV= 4.88 


Y = 55 


MAX-DMD= 130 


SIMULATION : 


AV E-COST= 37.89 


#0F 


SHOR ’ S= 


0. 1828 


SHSI MU= 6. 


73% 


SD-COST= 1.23 


SD- 


#OF SH'S 


= 0.0059 


NORMAL-APPROXIM. : 


A V E-COST= 36.71 


#OF 


SHOR’ S = 


0. 17 19 


SUMMARY STA 


TISTICS: 








SD/A VE= 


3. 25% 


PISIMA= 207.27$ 


PINORA= 


213.51$ 


DEVMON= 


3.12% 


D E7N SH= -5.96% 


DEVPI A= 


3. 01% 
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NO F DM D= 3 #OF TIME UNITS= 66 

MEAN 



REPLI=20000 I-S E ED= 394861 



#OF REP- 
CYCLES 



PROBA- 

BILITY 



UNIT-SHORT. 

LOST 



1: 


11 


1 : 


0.05 


1 


0. 


55 


1 


2: 


66 


2: 


0.30 


2 


19. 


80 


2 


3: 


33 


3: 


0. 10 


3 


3. 


30 


3 


>TAL 


M E AN = 


23. 


.65 


TOTAL 


S-DEV 


'= 4. 


17 



300.00 

300.00 

100.00 

Y= 30 



FIT 
N* P N*Q 

0.5 10.4 
19.8 46.2 
3. 3 29.7 

MAX-DMD= 110 



SIMULATION 
SHS IM U= 5.4615 

NORMAL— APPROXI M. : 



AV E— COST= 
SD— COST= 

A V E-COST= 



35.25 #OF S H O P. ' S = 
1.26 SD- #OF SH'S = 



31.08 #OF SHOE' S= 



0. 1292 
0.0046 

0. 1 142 



SUMMARY STATISTICS: 

SD/AVE= 3.57% PISIMA = 

DEVMON= -11.34% D EVNSH= 



272.73$ 
-1 1. 63% 



PINORA= 

DEVPIA= 



272. 09$ 
-0. 23% 



NOFDMD= 3 


#OF TIME 


UNI TS= 72 REPLI = 


20000 I-SS 


ED= 90 1941 


#OF REP- 


PROBA- 






UNIT-SHORT. 


FIT 


CYCLES 


BILITY 


MEAN 




COST 


N*P H*Q 


1: 18 


1: 0.20 


1: 


3.60 


1: 


300.00 


3.6 14.4 


2: 36 


2: 0.40 


2: 


14.40 


2: 


300.00 


14. 4 21.6 


3: 12 


3: 0.30 


3 : 


3.60 


3: 


100.00 


3.6 8.4 


TOTAL MEAN= 


21.60 


TOTAL S- 


DEV= 3.75 


Y= 23 


MA X-DMD= 66 


SIMULATION: AVE-COST = 


246.61 


#0 F SHOR'S = 


0.9046 


SHSIMU= 30. 


36% 


SD-COST= 


3.37 


SD 


- #OF SH'S 


= 0.0124 


NORMAL-APPROXIM. : AVE-COST = 


238.36 


#0 F SHOR'S = 


0. 8938 


SUMMARY STA 


.TISTICS: 












SD/A V E= 


1.36% 


PISIMA = 


272.63$ 




PINORA= 


266.67$ 


DEV MON= 


3.35% 


D EVNSH= 


-1.19% 




DEVPIA= 


-2. 19% 



NOFDMD= 3 #OF TIME UNITS= 40 REPLI=20000 I-5EED= 32467 



#OF REP- PROBA- 

CYCLES BILITY MEAN 




UNIT-SHORT. 

COST 


FIT 

N* P N*Q 


1: 10 1: 0.30 1: 
2: 20 2: 0.60 2: 
3: 40 3: 0.50 3: 


3.00 

12.00 

20.00 


1: 

2: 

3: 


100.00 

100.00 

300.00 


3.0 7.0 

12.0 8.0 
20.0 20.0 


TOTAL MEAN = 35.00 TOTAL S- 


DEV= 4. 


11 


Y= 38 


MA X-DMD= 70 


SIMULATION : AVE-COST = 
SHSIM U= 19.58% S D-COST= 


122. 37 

2.16 


# OF 
SD- 


SHOE' S= 
#0 F SH'S 


0.5349 
= 0. 0096 


NORMAL-APPROXIM. : AVE-COST = 


118. 98 


#0 F 


SHOE' S = 


0.5552 


SUMMARY STATISTICS: 

SD/A V E= 1.77% PISIMA= 

DEVHON= -2.77% DSVNSH= 


228. 78$ 
3.81% 


PINORA= 

DEVPIA= 


214. 29$ 
-6.34% 
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APPENDIX e 

SIHDLATION BESDLTS FOH FIVE DEMAND SOURCES 



NCFDMD= 5 


#OF TIME 


UNITS= 




60 REPLI = 


=20000 I-SE 


ED= 24 30 78 


#OF EEP- 


PROBA- 








UNIT-SHORT. 


FIT 


CYCLES 


BILITY 


MEAN 


COST 


N *P N*Q 


1: 10 


1; 0.10 


i 




1.00 


1 


200.00 


1.0 9.0 


2; 30 


2: 0.30 


2 




9.00 


2 


300.00 


9.0 21.0 


3: 20 


3: 0.50 


3 




10.00 


3 


100.00 


10.0 10.0 


4: 12 


4: 0.70 


4 




8.40 


4 


300.00 


8.4 3.6 


5; 15 


5: 0.90 


5 




13.50 


5 


200.00 


13.5 1.5 


TOTAL MEAN 


= 41.90 


TOTAL 


S 


- DE V= 4.01 




Y= 42 


MAX-DMD= 87 


SIMULATION: AVE-COST= 


333. 30 


#OF S HOR ’ S = 


1.5449 


SHS I M U= 43 


. 77% 


SD-COST= 


3.62 


SD-#OF SH'S 


= 0.0165 


NOBMAL-APPEOXIM. : AVZ-COST= 


336.42 


#0 F SHOE * S = 


1. 54 56 


SUMMARY STATISTICS: 














SD/A VE= 


1.09% 


PISIMA= 




215.74$ 




PINORA= 


217. 66$ 


DEVMON= 


0.93% 


DEVNSH; 




0.04% 




DEV P IA= 


0. 89% 




NOFDMD= 5 


ffOF TIME 


UNITS= 




60 R EPLI = 


=20000 I-SE 


ZD= 243078 


#OF REP- 


PROBA- 








UN IT- SHORI . 


FIT 


CYCLES 


BILITY 


MEAN 


COST 


N*P N*Q 


1: 10 


1: 0.10 


1 




1.00 


1 


200.00 


1.0 9.0 


2: 30 


2: 0.30 


2 




9.00 


2 


300.00 


9.0 2 1. 0 


3: 20 


3: 0.50 


3 




10.00 


3 


100.00 


10.0 10.0 


4: 12 


4: 0.70 


4 




8.40 


4 


300.00 


8.4 3.6 


5: 15 


5: 0.90 


5 




13.50 


5 


200.00 


13.5 1.5 


TOTAL MEAN 


= 41.90 


TOTAL 


S 


- DE V= 4.0 1 




Y= 46 


MAX-DMD= 87 


SIMULATION: AVE-COST= 


69.44 


#OF SHOR’S= 


0.3223 


SHS IM U= 12 


. 76% 


SD-CO ST= 


1.62 


SD— #OF SH’S 


= 0. 3074 


NORMAL-APPROXIM. : AVE-COST= 


69.07 


#OF SHOR* S = 


0.3173 


SUMMARY ST 


ATISTICS : 














SD/A V E= 


2.33% 


PISIMA; 




215. 13$ 




PINORA= 


217. 66$ 


DEVMON= 


-0.54% 


DEVNSH; 




-1. 70% 




DEVPIA= 


1. 13% 
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TT0FErMB~5 f OF - TXHZ - 'UNXTS=~6'U — — I = SFFU= — Z73075 - 



#0F EE 


?- PROBA- 


UNIT-SHORT. 


FIT 


CYCLES 


BILITY MEAN 


COST 


N*P N*Q 


1 : 10 


1: 0.10 1: 1.00 


1: 


200.00 


1.0 9.0 


2: 30 


2: 0.30 2: 9.00 


2: 


300.00 


9.0 21.0 


3: 20 


3: 0.50 3: 10.00 


3: 


100.00 


10.0 10.0 


4; 12 


4: 0.70 4: 8.40 


4: 


300.00 


8.4 3.5 


5: 15 


5; 0.90 5; 13.50 


5: 


200.00 


13.5 1.5 


TOTAL HE AN = 41.90 TOTAL S-DEV= 4.01 




Y= 48 


MAX-DMD= 87 


SIMULATION : AVI— COST= 24.09 


£OF SHOE ’ S= 


0.1130 


SHS IM U= 


5.08% SD-COST= 0.92 


SD 


- #OF SH’S 


= 0. 0042 


NORMAL- 


APPROXIM. : A V E-COST= 24.05 


#0 F SHOE' S = 


0. 1 105 


SUMMARY 


STATISTICS: 








SD/A V E= 


3.82% PISIMA= 213.08$ 




PINORA= 


217. 66$ 


DE VMO N= 


-0.15% DEVNSH= -2.25% 




DEVPIA= 


2. 15% 
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APPENDIX F 

OUTPUT TABLEAU OF THE SIMULATION PROGRAM FOR TWO DEMAND 

SOURCES 

+ ^^c*#^# + ****t***#**# **£:)[$***$* I*:):***#****#***#*#***#*##* *# # 

REPLIC' NS: 100300 1 TIME UNITS PEE PERIOD:202 I- SESD= 1 50449 3 

#OF REPCYCLS (1 ) : 4 4 P1=0.25 S MU1=1.00 6 PI 1= 100. 00 7 Y=7 8 

#0 F REPCYCLS (2) : 5® P2=0.75 10 MU2=3. 75 1 1 PI2=1000. 00 12 

RESULT OF BERNOULLI TRIALS FOR EACH REPAIR CYCLE: 

ESTIMATED MEAN OF SHORTAGE COSTS PER PERIOD: 13.43713 

ESTIMATED ST-DEV. FOR ESTIMATE OF THE MEAN: 0.3651® 

RESULT OF MONOVAR. APPROXIMATING FORMULA: 11.77115 

RESULT OF HYPERGEO. APPROXIMATING FORMULA: 9.3321® 

RESULT OF CDF APPROXIMATING FORMULA: 13.4391 7 

SUMMARY : i®MU= 4.751® SHEXA C= 1 . 36 1 2 o SHSIMU= 1.36%2i 

SD/MEAN=2.7%22 D EV MO N=- 12-4% DEVHYP=-30.6 % DEVCDF= 0 . 0 %2 3 

♦ ♦It****************************************************** 

Comments refe rr ing t o th e ou tput tableau: 

1. Number of simulated periods; 

2. Length of one period; 

3. Initial seed, necessary to reproduce identical simu- 
lation result ; 

4. Number of repair cycles for end item no. 1 per 
period. For instance, if the time period has a 
length of 20 time units and the number of repair 
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cycles is 4 then each repair cycle has a length of 5 
time units. 

5. Prolability of demand in one repair cycle for end 

item no. 1 ; 

6. dean demand for period for end item no. 1 ; 

7. Dnit shortage cost if one demand caused by end item 

no. 1 not satisfied; 

8. Purchased quantity available to satisfy period 
demand; 

9. Number of repair cycles for end item no. 2; 

10. Probability of demand in one repair cycle for end 

item no. 2 ; 

11. Mean demand for period for end item no. 2; 

12. Unit shortage cost, one demand caused by end item 
no. 2 not satisfied; 

13. Sample mean of the shortage costs per period as 
observed during simulation; 

14. Estimates precision of mean estimate. 

15. Monovariable (McMasters*) approximation formula; 

16. Hypergeometric approximation formula; 

17. CDF approximation formula; 

18. The last two lines provide a condensed simulation 
summary. 

19. Sum of expected demands from both demand sources for 
one period; 

20. -The program computes the distribution of the sum of 
the two binomial demand distributions and is there- 
fore able to output an exact expected fraction of 
periods in which shortages occur. In this example 
1360 periods with shortages ace the expected number 
so the expected fraction is 1.36%. 

21. Fraction of periods in which shortages occurred 
during simulation. When compared to SHEXAC it 
provides a numerical hint for the quality of the 
random numbers generated during the simulation. 
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22. The ratio of estimated standard deviation for the 
estimate of the mean to the sample mean. SD/MEAN 
decreases as the number of replications increases and 
the confidence interval for the mean derived from 
simulation becomes smaller. The value of SD/MEAN is 
the basis for the decision as to whether the number 
of replications ensures sufficient precision. The 
number of replications must be traded off against the 
desired precision. 

23. DEVMOK, DEVHYP, and DEVCDF express how far the three 

approximation formulas deviate from the simulation 
result. The complete last row allows a quick survey 
about the performance of the approximations. From 
the observed sample mean for the shortage costs per 
period and its estimated standard deviation, a confi- 
dence interval for the true unknown mean can be 
computed. It can be used to determine the value of 
the confidence level such that the confidence 

interval for the true mean would include the values 
of the approximation formulas. For instance if 
SD/MEAN = 5% and DEV MO N = -10%, it follows that the 
value of the monovariable approximation formula is 
two estimated standard deviations away from the 
sample mean and would be covered by a confidence 
interval with a confidence level' of at least -9545 
i.e. ?(|Z| < 2) =.9545 . 
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